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ABSTRACT 

The low-mass X-ray binary 4U 0614+091 is a source of sporadic thermonuclear (type I) X-ray bursts. We find bursts with a wide 
variety of characteristics in serendipitous wide-field X-ray observations by the WATCH on EURECA, the ASM on RXTE, the WFCs 
on BeppoSAX, the FREGATE on HETE-2, the IBIS/ISGRI on INTEGRAL, and the BAT on Swift, as well as pointed observations with 
the PCA and HEXTE on RXTE. Most of the bursts are bright, i.e., they reach a peak flux of about 15 Crab, but a few are weak and 
only reach a peak flux below a Crab. One of the bursts shows a very strong photospheric radius-expansion phase. This allows us to 
evaluate the distance to the source, which we estimate to be 3.2 kpc. The burst durations vary generally from about 10 sec to 5min. 
However, after one of the intermediate-duration bursts, a faint tail is seen to at least about 2.4 hours after the start of the burst. One 
very long burst was observed, which lasted for several hours. This superburst candidate was followed by a normal type-I burst only 
19 days later This is, to our knowledge, the shortest burst-quench time among the superbursters. The observation of a superburst in 
this system is difficult to reconcile if the system is accreting at about 1% of the Eddington limit. We describe the burst properties in 
relation to the persistent emission. No strong correlations are apparent, except that the intermediate-duration bursts occurred when 
4U 0614+091 's persistent emission was lowest and calm, and when bursts were infrequent (on average roughly one every month to 
3 months). The average burst rate increased significantly after this period. The maximum average burst recurrence rate is about once 
every week to 2 weeks. The burst behaviour may be partly understood if there is at least an appreciable amount of helium present in 
the accreted material from the donor star. If the system is an ultra-compact X-ray binary with a CO white-dwarf donor, as has been 
suggested, this is unexpected. If the bursts are powered by helium, we find that the energy production per accumulated mass is about 
2.5 times less than expected for pure helium matter. 

Key words. Accretion, accretion disks - binaries: close - Stars: individual: 4U 0614+091 - Stars: neutron - X-rays: binaries - X-rays: 
bursts 



1. Introduction 

Type I X-ray bursts (Grindlay et al. 1975, Belian et al. 1976, 
Hoffman et al. 1978; hereafter bursts) result from thermonuclear 
shell flashes on a neutron star, which is caused by the ignition 
of either He and/or H-rich material supplied by a binary com- 
panion star (Hansen & van Horn 1975, Woosley & Taam 1976, 
Maraschi & Cavaliere 1977, Lamb & Lamb 1978; for reviews 
see Lewin et al. 1993, Strohmayer & Bildsten 2006). Bursts gen- 
erally appear as short transient events wherein the X-ray inten- 
sity rises rapidly on a time scale of seconds, and decays in an 
exponential fashion back to the pre-burst level. The decay lasts 
almost always longer than the rise. Burst durations range from 
several seconds up to half an hour. The burst spectra generally 
harden during the rise and soften during the decay. This has been 



attributed to the heating and cooling of the uppermost layers of 
the neutron star. The spectra can be satisfactorily described by 
black-body emission from spherical regions with radii of around 
10km at inferred temperatures up to ^r^3keV. The burst-to- 
burst time intervals are typically of the order of hours to days. 
They are thought to be determined by the time for the neutron 
star to accumulate enough fuel to power another burst. 

Under certain conditions, the local luminosity may reach or 
exceed the Eddington limit and matter may be pushed outward. 
As a result, the neutron star photosphere also moves outward. 
Consequently, the emitting area increases and the observed in- 
ferred black-body temperature drops. When the surge of energy 
release is over, the photosphere gradually returns to its pre-burst 
radius. During this phase, the emitting area decreases and the 
infeiTed temperature increases. During the expansion and con- 
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traction phase, the luminosity is expected to be close to the 
Eddington limit. After the photosphere recedes to its pre-burst 
radius (called 'touch-down') cooling is typically observed. Such 
bursts are referred to as photospheric radius-expansion type I X- 
ray bursts (e.g., Lewin et al. 1984, Tawara et al. 1984; hereafter 
radius-expansion bursts). 

A few bursts last hours to half a day (Cornelisse et al. 2000, 
Strohmayer & Brown 2002; for a review see, e.g., Kuulkers 
2004). Such bursts show the same characteristics as the typ- 
ical bursts described above, except that the energy output is 
about a thousand times higher. They recur on time scales of 
months to years. Their origin is different from the typical bursts, 
because they are thought to be due to unstable burning of C 
deeper in the neutron star envelope (Gumming & Bildsten 2001, 
Strohmayer & Brown 2002). These long events are called super- 
bursts (Wijnands 2001). Typical bursts are seen up to the time 
of the superbursts. However, they cease to occur after the super- 
burst, for up to about a month (see, e.g., Cornelisse et al. 2002a, 
Kuulkers et al. 2002a). This is thought to be due to the super- 
burst heating up the surface layer, and therefore preventing the 
unstable ignition of H and/or He (Gumming & Macbeth 2004). 

There are also bursts that have durations and energy releases 
intermediate between typical bursts and superbursts. They are 
referred to as intermediate-duration X-ray bursts (Gumming et 
al. 2006; hereafter intermediate-duration bursts). They show de- 
cay times ranging from several to tens of minutes and energy 
outputs of about 10"*' erg. Most of these events occur in sources 
with low persistent luminosity (Lpers^O.OlLEdd, where LEdd rep- 
resents the Eddington luminosity) and are thought to be due to 
flashes of relatively thick He layers in ultra-compact X-ray bi- 
naries (UGXBs; i.e., low-mass X-ray binaries with an orbital pe- 
riod less than an hour), see in 't Zand et al. (2005), Gumming 
et al. (2006), Falanga et al. (2008), and references therein. The 
intermediate-duration bursts in the H-rich and luminous system 
GX 17-1-2 (Kuulkers et al. 2002b) probably have a different ori- 
gin. 

The low-mass X-ray binary 4U 0614-1-091 has long been 
known to be a source of bursts, but this identification was based 
on the detections of only several of them, with relatively poor, 
position determination^ as well as coarse timing and/or X- 
ray spectral information (Lewin 1976, Swank et al. 1978, Brandt 
et al. 1992, 1993a,b, Brandt 1994, Brandt & Lund 1995). More 
recently, an hours-long flare was reported, which showed charac- 
teristics resembling superbursts (Kuulkers 2005; see Sect. 13.21 ). 
Brandt et al. (1992) used the observation with GRANAT/WATCH 
of a peak burst flux of ~ 10 Grab (6-20 keV) to infer, assuming 
this peak flux does not exceed the Eddington limit for a 1 .4 M© 
neutron star with a He-rich atmosphere, that the source distance 
must be about 3.2kpc. 4U 0614+091 is, therefore, one of the 
closest to Earth of the known actively bursting systems (see, e.g., 
Jonker & Nelemans 2004, Galloway et al. 2008). 

In the last decade, HETE-2 triggered thirteen times on 
bursts from 4U 06 14+091. Moreover, three strong bursts from 
4U 0614+091 triggered the y-ray burst monitors on board 
INTEGRAL (see Ghelovekov et al. 2007) and Swift. Strohmayer 
et al. (2008) discovered burst oscillations near 415 Hz in one of 
the brightest ^vvZ/f/BAT bursts. 

It is difficult to understand 4U 06 14+091 as the origin of 
bursts. Since it is thought to be an UGXB (in 't Zand et al. 2007, 
Shahbaz et al. 2008, and references therein), it is supposed to 
be so compact that the donor star can only be a non-degenerate 



' Note that spatially the closest known X-ray burster to 4U 0614+091 
is MXB 05 1 3-40, with a distance of 5 1° . 



H-deficient star or a white dwarf. Evolutionary models suggest 
that the accreted material could be rich in G/O or He (Nelemans 
et al. 2009). However, no evidence of He (or H) has been found 
so far in optical spectra, with rather stringent upper limits: He 
and/or H are at most present at the 10% level (Werner et al. 
2006). These spectra suggest the donor star to be a G/O white 
dwarf (Nelemans et al. 2003, 2006; Werner et al. 2006). This 
seems to be at odds (see Juett et al. 2001, Juett & Ghakrabarty 
2003) with the fact that we see bursts in 4U 0614+091, which 
are usually understood as being due to unstable ignition of He 
or H. Moreover, 4U 0614+091 is an interesting case among the 
sources showing superbursts. Superbursts are expected to occur 
only in systems with accretion rates in excess of about 10% 
of the Eddington rate, i.e., the rate which yields emission at 
the Eddington limit (Gumming & Bildsten 2001, Strohmayer & 
Brown 2002). However, the low persistent X-ray luminosity of 
4U 06 14+091 indicates an accretion rate which is an order of 
magnitude lower (e.g.. Ford et al. 2000). 

In this paper we report on the detection of bursts from 
4U 0614+091 made with various instruments. In Sects. ITTI and 
12.21 we describe the instruments and the data analyses of the 
bursts seen between 1992 and 2007, as well as the analysis of the 
persistent emission. We complement this with information ex- 
tracted for the bursts seen before 1992 (Sect. 123] ). An overview 
of all the typical and intermediate-duration bursts seen appears in 
Sect. 13.11 We collectively refer to both types of bursts as normal 
bursts. One of them is a radius-expansion burst; it is discussed 
in Sect. 13.1.21 We then describe the results of our analysis of the 
superburst (Sect. 13.21 ). our burst oscillation search (Sect. 13.3b . 
and the characteristics of the persistent emission of the source 
(Sect. I3.4l i. Finally, we show the long-term X-ray behaviour of 
4U 0614+091 (Sect. [331 ) and close the paper with a discussion 
of our results (Sect.|4|. 

2. Observations and data analysis 

Almost all the observations between 1992 and 2007 described 
in this paper were serendipitously obtained using six different 
X-ray instruments on just as many space-borne observatories, 
see Table [T] Only the observations with the RXTE/PCA and 
RXTE/HEXTE were purposely dedicated to 4U 06 14+091. In 
the next subsections we describe these instruments and the anal- 
ysis of the data in more detail. We usually refer to the instru- 
ments by acronyms only. In Table [T] we also give the total time 
of exposure for observations with 4U 0614+091 in the field of 
view and the total number of bursts detected from this source in 
these observations. 

2.1. Observations and instrument-specific data analysis 
issues 

2.1.1. EURECA\NAJC\-\ 

EURECA (European REtrievable CArrier) carried the Wide 
Angle Telescope for Cosmic Hard X-rays (WATGH; Lund 1985, 
Brandt et al. 1990). WATGH operated from August 1992 to June 
1993. It comprised a rotating modulation collimator (RMG) and 
yielded images through a cross-correlation method. It viewed a 
circular field with a radius of about 65°, was sensitive between 
6 and 150keV, and had an effective area of about 45 cm^. For 
events with a duration longer than the rotation period of the 
RMG (^1 s), the relative position accuracy was better than 1°, 
and, in favourable cases, was as good as about 0.1°. The flux 
sensitivity was about lOOmGrab in one day. 
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Table 1. Overview of instruments that have detected bursts from 
4U 0614+091 between 1992 and 2007. 

Time span Satellite/Instrument Escm" ^exp^ 

(keV) (days) 

Aug 1992- Jun 1993 EURECA/WATCH 6-150 -50 3 

Dec 1995 - Aug 2007 RXTE/ASM 1.5-12 ^58 7 

Dec 1995 - Aug 2007 RXTE/PCA 2-60 ^23 2 

Dec 1995 - Aug 2007 RXTE/HEXTE 15-250 ==23 2 

Apr 1996- Apr 2002 BeppoSAX/W¥C 2-28 ==27 1 

Oct 2000 - Mai- 2007 //£r£-2/FREGATE 6^00 -194 13 

Oct 2002 - Aug 2007 INTEGRAL/lSGRl 15-1000 ^22 2 

Oct 2002 - Aug 2007 INTEGRAL/mM-X 3-35 =0.8 

Nov 2004 - Aug 2007 Swift/BAT 15-150 =26'' 2 

" Nominal sensitive energy range. 
* Total observation net exposure time. 

Number of bursts seen. 
'' Total effective exposure time (corrected for off-axis response). 



analysis of the coded-aperture data and is made available via the 
internet, on a dwell by dwell basis as well as 1-day averagesQ 

The intensity history of 4U 06 14+091 up to September 2007 
was obtained from each of 62720 dwells, more or less evenly 
spread over the years (see also Sect. B.S.TT l. We used the follow- 
ing method to search for bursts: 1) identify candidate dwells in 
which the source intensity minus its uncertainty was larger than 
the overall mean intensity of 4U 0614+091 by 4 times the rms 
uncertainty in the overall mean intensity; 2) visually inspect raw 
count rates at 1-s time resolution within the candidate dwells for 
evidence of burst-like temporal behaviour; 3) use the hardness 
ratios, defined as the time-bin-wise ratios of the count rates in 
the 5-12 keV band to those in the 1.5-5 keV band, to verify the 
typical hard-to-soft behaviour of bursts. Note that this technique 
is not suitable for bursts that last substantially longer than a dwell 
or for bursts that do not start within a dwell. Our search yielded 
6 bursts and one superburst. 



Aug 1992 -Jun 1993 
Dec 1995 - Aug 2007 
Dec 1995 - Aug 2007 
Dec 1995 - Aug 2007 
Apr 1996- Apr 2002 
Oct 2000 - Mar 2007 
Oct 2002 - Aug 2007 
Oct 2002 - Aug 2007 
Nov 2004 - Aug 2007 



EURECA/WATCH 
RXTE/ASM 
RXTE/PCA 
RXTE/HEXTE 
BeppoSAX/W¥C 
HETE-2/FREGATE 
INTEGRAL/lSGRl 
INTEGRAL/mM-X 
Swift/BAT 



6-150 


-50 


1.5-12 


=58 


2-60 


=23 


15-250 


=23 


2-28 


=27 


6^00 


-194 


15-1000 


=22 


3-35 


=0.8 


15-150 


=26'' 



4U 0614+091 was located in the ^1 steradian field of view 
from mid January 1993 until the end of April 1993, with an ob- 
serving efficiency of a little more than 40%. EURECA/WATCH 
observed a total of three bursts during that time. Preliminary re- 
ports of these events can be found in Brandt et al. (1993a,b), 
Brandt (1994) and Brandt & Lund (1995). 

For our analysis, we used count rates as a function of time 
along with vignetting estimates for the source. Background rates 
were determined by polynomial fits to the count rates outside the 
intervals of each burst. We define the 3cr confidence level above 
which we regard a signal as significant as 3 times the square root 
of the number of background counts in each time bin. 

2.1.2. RXTEAW-Sky Monitor 

2.1 .2.1 Scanning Sliadow Cameras. The Rossi X-ray Timing 
Explorer (RXTE; launched December 1995) All-Sky Monitor 
(ASM) consists of three Scanning Shadow Cameras (SSCs, 
hereafter called SSCO, SSCl and SSC2; Levine et al. 1996). 
Each SSC views a 12°xll0° (full-width at zero response, 
FWZR) field through a random-slit coded mask. The field of 
view of one SSC is 90° from the colocated field centres of the 
other two SSCs. The assembly holding the three SSCs is gener- 
ally held stationary for a 90s "dwell". A drive then rotates this 
assembly through 6° between dwells. This yields good sky cov- 
erage: as much as 80% of the X-ray sky is covered every 90min 
orbit around the Earth. 

For each dwell, histograms of counts ('raw' counts, i.e., not 
corrected for transmission, sources in the field of view, back- 
ground, etc.) as a function of position in each detector are 
recorded for three energy bands (roughly 1 .5-3, 3-5, 5-12 keV). 
Also, the total number of counts registered in each SSC are 
recorded in 1/8 s time-series bins in the same energy bands. For 
observations done before March 2001 (MJD 51970), there is 
imaging data integrated over entire dwells, and count rate data 
with no imaging information. Since March 2001 event-by-event 
data are telemetered so that both temporal and imaging informa- 
tion can be extracted for portions of dwells. For sources away 
from bright sources, the sensitivity is roughly 10-15 mCrab (lo- 
in the 1.5-12 keV band) for a reliable single-dwell source detec- 
tion; multi-day averaging improves the sensitivity down to about 
a few mCrab. The average intensity of each known source in the 
field of view is determined by the instrument team through an 



2.1.2.2 (Time-resolved) X-ray spectral analysis of ASM data. 
The 3-channel ASM data can constrain simple spectral models 
like power laws and black-body radiation (see, e.g.. Ford et al. 
1996, Kuulkers 2002, Keek et al. 2008). In our spectral analyses 
we employed the same method as in Keek et al. (2008). This 
includes calibrations of the source flux with average ASM data 
of the Crab over 200-day intervals around each data point. 

The spectral analysis can also be satisfactorily applied to raw 
ASM data when dealing with bursts. If one assumes that the 
burst emission is not influenced by persistent emission (see, e.g., 
Kuulkers et al. 2002b), net-burst count rates can be obtained by 
subtracting the average raw pre-burst count rates from the raw 
burst count rates. If one also knows the celestial position of the 
burst, the effective mask/collimator transmission can be com- 
puted. This allows the estimation of equivalent on-axis net burst 
count rates for the given energy bands and these, after renor- 
malization with the Crab as described above, can be fit with a 
black-body spectrum model. Whenever a burst was simultane- 
ously observed by two SSCs we combined the renormalized net 
count rate information from the two cameras. We derived the 
three-channel transmission-corrected net-burst spectra from the 
1-s raw count rates. We calculated the errors by taking square 
roots of the numbers of counts per each time bin. We started 
with a spectral time bin size of 1 s (2S 0918-549, see below) or 
2 s (4U 0614+091). Each time the net-burst rates foUowing the 
peak decreased by a factor of V2, we doubled the spectral time- 
bin size. 

Since all the corrections may introduce errors in our analysis 
of the raw ASM data, we first applied our method to bursts from 
2S 0918-549 that, like 4U 0614+091, lies in a relatively empty 
field in the sky (spatially the closest bursters to 2S 0918-549 are 
EXO 0748-676 at ^16° away and the transient GS 0836-429 at 
^14° away). 2S 0918-549 is similar to 4U 0614+091 in that it is 
also thought to be an UCXB which shows low persistent accre- 
tion and infrequent bursts (see, e.g., in 't Zand et al. 2005). Most 
of the 2S 0918-549 bursts show a radius-expansion phase. By 
comparing the ASM time-resolved X-ray burst spectral analysis 
to that derived from high-quality X-ray spectra observed with 
well-calibrated detectors, one can verify our ASM burst-spectral 
analysis method. 

We searched the ASM data on 2S 0918-549 for bursts us- 
ing the same algorithm as that used for 4U 06 14+091. We found 



^ ihttp : //heasarc . gs fc.nasa.gov/docs/xte/asm.products.htinl| 
and ltrttpTZ/xteTiTt .edu/ASM J.C ■ htinl| 
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5 bursts (see Table O over the more than ^11 year time span 
of observations, with a net effective exposure time of 47.8 days 
between 1996 January and 2007 August. Three of them were al- 
ready noted by in 't Zand et al. (2005); our third X-ray burst, 
which occurred on UT 2002 August 23, was not noted by them 
probably due of its weakness (because of relatively low trans- 
mission, see Table [3]l. By comparing our ASM time-resolved 
spectral fit results with those for the bursts seen with the WFC 
(see in 't Zand et al. 2005) and PCA (see Galloway et al. 2008), 
we find comparable values for all spectral parameters. The ASM 
spectral fits show similar kind of radius expansion phases which 
reach similar derived peak fluxes as the other instruments. We, 
therefore, conclude that our time resolved X-ray spectral analy- 
sis of the raw ASM data gives consistent results with that seen 
from other instruments, and therefore our method can be trusted. 

2.1 .3. RXTE PCA and HEXTE 

The Proportional Counter Array (PCA; Bradt et al. 1993, Jahoda 
et al. 2006) onboard RXTE provides a large collecting area (max- 
imum net geometric area of about 8000 cm^) and high time 
resolution (down to fis). It consists of 5 proportional counter 
units (PCUs) behind 1° (full-width at half maximum, FWHM) 
collimators. It is sensitive in the 2-60 keV range down to 
^0.2mCrab. 

The High-Energy X-ray Timing Experiment (1 5-250 keV; 
Rothschild et al. 1998) onboard the same satellite consists of 8 
detectors with a total area of about 1600cm^. The 8 detectors 
are split up in two clusters of 4 each. During normal operations, 
each cluster is alternately pointed on and off the source generally 
every 16 or 32 s, to provide near real-time background measure- 
ments (note that at the end of 2006 cluster A was fixed to always 
view the source). HEXTE can measure a typical X-ray source 
down (at the 3cr level) to about 1 mCrab up to ^lOOkeV in 10^ s. 

We first inspected the light curves at 1 s time resolution of 
all the publicly available PCA data up to September 2007 on 
4U 06 14-1-091 that was collected in the Standard 1 mode. A cou- 
ple of strong burst-like events were seen; all except one were 
due to detector-related events. The one celestial event occurred 
simultaneously with the burst seen by FREGATE, on UT 2001 
February 4 (MJD 51944; see Sect. l2X5T l. Only the first 30 s 
of the event were covered, after which the PCA was shut off 
because the count rate (>55kcts ') exceeded the 'High Rate 
Monitor' safety threshold. Apart from the first 30 seconds of the 
burst, the PCA also observed the tail of the burst. In observations 
taken on UT 2000 September 2 (MJD 5 1789) we found the tail 
of another burst for which the onset was missed. 

We used version 6.5 of the HEASOFT software suite for 
our RXTE data analysis. For the time-resolved analysis of the 
prompt burst on MJD 5 1944 we used the PCA data from either 
Event mode or Burst Catcher mode when available (the high 
count rates reached during the maximum of the burst resulted 
in data losses near the ends of the 1 s buffers of the Event mode). 
These modes provided numbers of counts in 64 channels cover- 
ing the PCA energy range at time resolutions of 122 fj.s and 8 ms, 
respectively. At the time the prompt burst occurred, 4 PCUs 
(PCU 0-3) were on; the data from all layers and all PCUs were 
automatically combined in these modes. During the spike (see 
Sect. 13.1.2b we created time-resolved spectra at 2 ms resolution. 
During the rest of the prompt burst, we used a time resolution of 
0.125 s. Although the majority of the prompt burst was observed 
by HEXTE, we do not use the HEXTE data for the burst spectral 
analysis because data from a broader energy range (7-400 keV) 
was obtained by FREGATE (see Sect. 12.1.5b . For the tail of the 



bursts on MJD 51789 and MJD 51944 we used the PCA data 
from the Standard 2 mode; this mode provides counts in 129 
channels covering the PCA energy range with a time resolution 
of 16 s. We formed spectra averaged over 112 s intervals from the 
events occurring in all layers of the PCUs which were on during 
and before or after the burst; this was PCU 3 on MJD 51944 and 
PCUs 2,3 on MJD 51789. 

Dead-time correction is only possible for all spectra with a 
time resolution equal or higher than 0.125 s. We subtracted the 
pre-burst persistent emission from the burst emission. Our time- 
resolved burst spectral fits of PCA data covered the 3 to 20 keV 
energy band. We included a 1 % systematic uncertainty in each 
spectral bin in addition to the usual statistical uncertainties. 

We extracted information on the persistent emission only 
when RXTE data was available within 1 day of a burst (see 
Sect. [IS- We did this for data from both the PCA and HEXTE 
by applying standard criteria, i.e., by filtering out data taken 
at elevations less than 10° or with an offset from the source 
greater than 0.02°. Also, we created instrument response and 
background files following the standard analysis threads using 
the latest information available^ We corrected the count-rate 
spectra for dead-time losses. For the PCA we only used PCU 2 
data (across the RXTE mission this PCU has the highest duty 
cycle) obtained with the Standard 2 mode. Pulse-height spectra 
were formed from the events from all layers and we included a 
1% systematic uncertainty in each spectral bin. We performed 
fits to the persistent emission in the 3-30keV and 17-100 keV 
bands for the PCA and HEXTE, respectively. A multiplicative 
constant was applied to the HEXTE spectra to account for the 
uncertainty in the relative normalizations of the instruments, i.e., 
both HEXTE Cluster A and B were allowed to vary with respect 
to the PCA. This multiplicative factor was found to be between 
0.6 and 0.8. 

For the comparison (see Sect. 13.1.2b of the burst observed 
on MJD 51944 with that of the superbursts of 4U 1820-303 
(Strohmayer & Brown 2002) and 4U 1636-536 (Strohmayer 
& Markwardt 2002) at hard X-ray energies, we extracted light 
curves from data obtained by HEXTE in the 15-60keV band. 
In this process we applied the same standard filtering crite- 
ria as those described above. We corrected the light curves for 
dead time and background following the standard procedure^ 
and normalized them to the Crab count rate in the same en- 
ergy band. For this analysis we only used the Standard Modes 
(Archive Spectral Bin, 64-bin spectra produced every 16 s) data 
from Cluster A. 

To characterise the overall spectral behaviour of 
4U 0614+091 we created so-called colour-colour and hardness- 
intensity diagrams (CD and HID, respectively). All publicly 
available pointed PCA observations of 4U 0614+091 up to 
September 2007 were used. We extracted background and 
dead-time corrected count rates from the Standard 2 data using 
the same standard procedures and filtering criteria as applied 
to the spectral data, with 16 s time resolution. We used the 
following energy bands: 2.0-3.5 keV (A), 3.5-6.0keV (B), 
6.0-9.7 keV (C), and 9.7-16.0keV (D). We define the soft and 
hard colours (SC, HC) as the ratios of the count rates in the 
various bands: SC=B/A and HC=D/C; the intensity (Int) is 
defined as the sum of the count rates in the four energy bands: 
Int=A+B+C+D. We normalized both the colours and intensity 
to the Crab values nearest in time to account for changes in 
detector gas gain (see, e.g., Kuulkers et al. 1994, van Straaten 



^ http : //heasarc . nas a . gov/ docs/xte/data^nalysis . html| 
|http : //heasarc . gsf c . nasa . gov/docs/xte/recipes/hexte . html| 
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et al. 2003). Since the values do not change significantly within 
an observation (i.e., on time scales of order an hour), they were 
averaged over one observation. 

2.1.4. Bep/7o5AX Wide-Field Cameras 

The Wide-Field Cameras (WFCs; Jager et al. 1997) were two 
identical coded-aperture instruments onboard the BeppoSAX 
satellite (Satellite per Astronomia X; Boella et al. 1997), which 
were operated between April 1996 and April 2002. The field of 
view was 40°x40° (FWZR), the angular- resolution 5' (FWHM) 
and the source-location accuracy was generally better than 1' 
(99% confidence). The detectors were sensitive to the energy 
range 2 to 28keV and had a net collecting area of 140cm-. The 
on-axis detection threshold was of the order of 0.3 Crab for a 1 s 
observation and a few mCrab for a 10^ s observation. 

The WFCs pointed in opposite directions with respect to 
each other and perpendicular to the Narrow-Field Instruments 
(NFIs) on the same satellite. Since the pointing directions of 
the WFCs were usually governed by the observations of the 
NFIs, the WFC sky coverage was not uniform. 4U 0614-1-091 
was mostly seen at large ofF-axis angles. 

We extracted the reconstructed source flux in the 2-25 keV 
bandpass with a time resolution of 2 s whenever 4U 0614-)-091 
was in the field of view. We also extracted information on all 
photons detected on the appropriate detector with varying time 
resolutions between 0.5 and 8 s. Light curves were formed from 
both the fluxes and the photon rate data. Bursts were searched 
for by eye as well using an automatic algorithm (see Cornelisse 
et al. 2003). We found only one burst. 

For the generation of WFC X-ray spectra we first cross- 
correlated the detector data with the expected imaging response 
of WFC unit 2 for 4U 0614+091 (see in 't Zand 1992 for more 
detail). Background radiation is then automatically subtracted. 
We then extracted 2-28 ke V spectra from the resulting imaging 
data. Burst time bins were defined such that the significance of 
the source (i.e., the photon flux divided by the standard deviation 
expected from all other photon sources) is at least 10. Experience 
shows that spectra are ill defined for lower significances. The 
X-ray spectrum of the persistent emission around the burst was 
determined from the whole observation in which the burst oc- 
curred (with an exposure time of 14.3 ksec). This spectrum was 
subtracted from the burst emission in our time-resolved burst 
spectral analysis. The burst and persistent spectra were fit in the 
2-28 keV band. 

2.1.5. HETE-2 FREGATE 

One of the instruments onboard the High Energy Transient 
Explorer satellite {HETE-2; Ricker et al. 2003) that operated be- 
tween 2000 and 2007, was the omnidirectional y-ray spectrom- 
eter named French Gamma Telescope (FREGATE; Atteia et al. 
2003). FREGATE consisted of 4 detectors, mounted in two pairs. 
It was sensitive in the 6-400 keV range, with a maximum effec- 
tive area of 158 cm^. The field of view was 70° (half- width at 
zero response, HWZR), but it had no imaging capacities. It pro- 
vided continuous 128 channel energy spectra with 5 s time reso- 
lution (except in the beginning of the mission when it was 10 s), 
as well as continuous four-channel energy spectra with 0.164 s 
time resolution (0.327 s in the beginning). 

A burst trigger occurred when there were two coincident 6cr 
or higher excesses on two of the FREGATE detectors within a 
time bin; this corresponds to an increase by 1.2 Crab or higher 



in 5.24 sec for a burst observed on-axis in a source-free region. 
When this happened, 256000 individual photons were time- and 
energy-tagged. FREGATE triggered 13 times on bursts coming 
from the direction of 4U 0614+091 (see Table[T]l. Three of them 
were weak and exhibited low signal to noise (on MJD 52961, 
MJD 53074 and MJD 53740). The first two FREGATE events 
(MJD 5 1944 and MJD 52322) were already reported by Barraud 
(2002); they occurred when the Wide Field X-ray Monitor 
(WXM, imager onboard HETE-2, 2-25 keV; Shirasaki et al. 
2003) was not operating. The remaining bursts were detected 
by the WXM, and some of these were also detected by the Soft 
X-ray Camera (SXC, CCD-based imager onboard HETE-2, 0.5- 
2 keV; Villasenor et al. 2003). The WXM and SXC, when operat- 
ing, were able to get precise burst localizations (see, e.g., Suzuki 
et al. 2004). Two additional bursts were seen by the WXM, 
but FREGATE was not operating at those times (MJD 53041 
and MJD 54101). Due to the anti-solar pointing of HETE-2, 
4U 0614+091 was within the field-of-view of FREGATE dur- 
ing nearly 4 months every year. We estimate the total exposure 
on 4U 0614+091 from 2001 to 2006 to be about 194 days. 

The extraction of burst spectra relies on a quadratic fit of 
the 'background' (i.e., true background plus any emission from 
sources in the field of view) during tens of seconds before and af- 
ter a burst in all the 128 energy channels. The availability of this 
simple background spectrum allows for the subtraction of the 
background channel by channel. We used the counts above the 
fitted background to build the count spectrum of the burst, which 
we then deconvolved using the response matrix constructed with 
the known gain and angular response of the detector We per- 
formed the X-ray spectral fits to the integrated and time-resolved 
(10s resolution for the burst on MJD 5 1 944, 5 s for the rest) burst 
emission in the well-calibrated 7-40 keV range. For several of 
the bursts we also extracted spectral information at a higher time 
resolution in a similar way. 

2.1.6. INTEGRAL IBIS/ISGRI and JEM-X 

One of the two main instruments onboard INTEGRAL 
(International Gamma-Ray Astrophysics Laboratory; Winkler 
et al. 2003; launched October 2002) is IBIS (Imager on Board 
the INTEGRAL Satellite; Ubertini et al. 2003). It comprises two 
detector planes that view the sky through a coded mask. The field 
of view is 29°x29° (FWZR) and the angular resolution is 12' 
(FWHM). We use data collected with one of its detectors: the 
INTEGRAL Soft Gamma-Ray Imager (ISGRI) which is sensi- 
tive in the ^15 keV to 1 MeV range with a total effective area of 
about 2600 cm^ (Lebrun et al. 2003). 

One of the two supplementing monitors onboard INTEGRAL 
is the Joint European X-ray Monitor (JEM-X; Lund et al. 2003). 
JEM-X consists of 2 identical units, which are both sensitive in 
the 3-35 keV band. Most of the time only one unit is operat- 
ing. The angular resolution is 3' (FWHM); one single unit has 
a detector area of about 500 cm^. The units have a circular view 
with diameter of about 13° (FWZR), i.e., narrower than ISGRI. 
In practice, the transmission of the collimator beyond an ofF-axis 
angle of 5° is so low that only the brightest sources can be ob- 
served at larger angles. JEM-X is sensitive to X-ray bursts seen 
on-axis down to about 0. 1 Crab at the 5cr level in 5 s. 

The INTEGRAL Burst Alert System (IB AS; Mereghetti et 
al. 200313 is the automatic on-ground software that in near-real 
time searches for y-ray bursts in IBIS data and promptly publicly 
distributes results. For a trigger time interval of 1 s and a certain 

^ [http : //ibas . iasf-milano . inaf . it/| 
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threshold value (currently about 8cr), a minimum flux of about 
0.5-0.75 Crab (20-200 keV) is required to trigger a typical j- 
ray burst (and to produce enough counts to locate the position in 
the deconvolved image). One burst from 4U 0614+091 triggered 
IBAS on UT 2005 March 31 (nr 2441 on MJD 53460). 

To search offline for bursts and to study the hard X-ray long- 
term persistent behaviour of 4U 0614+091, we analysed all pub- 
lic INTEGRAL pointing data in the ISOC Science Data ArchiveQ 
Up to September 2007 there were ^840 pointings available 
where 4U 06 14+091 was within the ISGRI field of view with a 
total exposure time of about 25 days. Since most of these ob- 
servations were taken with the prime objective of instrument 
calibration on the Crab, many pointings (about 200) were car- 
ried out in non-standard modes. The total exposure time of all 
ISGRI standard pointings is about 22 days. All public JEM-X 
data whenever 4U 0614+091 was in the field of view (i.e., <6° 
ofi'-axis) amounts to a total exposure of only about 70ksec. In 
our analysis we use the high-energy source ISDC reference cat- 
alog (see Ebisawa et al. 2003; v. 27) as input source catalog. 
We processed the ISGRI and JEM-X data using the latest avail- 
able Off-line Scientific Analysis software (OSA; see Courvoisier 
et al. 2003), v. 7.0. The description of the algorithms used in the 
ISGRI and JEM-X scientific analyses can be found in Goldwurm 
et al. (2003) and Westergaard et al. (2003), respectively. 

We I) processed the ISGRI 15-20keV band data using 
OSA with default parameters, 2) analysed the data through to 
the imaging step, and 3) extracted light curves for all detected 
sources with 10 s time resolution^ Next, a potential onset of a 
burst was flagged when, in a time bin, the difference between the 
source count rate and the average source count rate in the whole 
pointing exceeded four times the standard deviation of the count 
rates in the whole pointing. The count rate versus time was then 
examined around each flagged time to check for the presence of a 
shape consistent with that of a type I X-ray burst. If this was the 
case, we generated reconstructed images within the good-time 
interval covering the whole burst, and checked them visually to 
verify that the event indeed originated from 4U 06 14+09 1 . 

Using the above described procedure we confirmed the 
burst which triggered IBAS in 2005 (see above). We also de- 
tected another burst from 4U 06 14+091 on UT 2003 August 16 
(MJD 52867). Unfortunately, the latter burst was seen far off- 
axis, about 16.5°, near the edge of the detector (see Fig. IA.3b . 
and thus again outside the field of view of JEM-X. We found no 
bursts from 4U 06 14+091 in the JEM-X data. 

For the time-resolved spectral analysis of the 2005 burst we 
divided the burst only in three parts so that the rise, the top and 
the decay were covered with satisfactorily spectral quality, in 5 
spectral energy bins covering 17 to 40keV. The statistical qual- 
ity of the 2003 burst precludes a time-resolved analysis. The 
ISGRI spectrum of the persistent emission around the 2005 burst 
is derived from the 1800 s single pointing during which the burst 
occurred. The net-burst spectra were calculated by subtracting 
the persistent spectrum from the burst spectra. X-ray spectral 
fits to the burst and persistent ISGRI spectra were performed in 
the 17-40 keV and 20-200 keV bands, respectively. During the 
INTEGRAL observations around the 2003 burst 4U 06 14+091 
was always far off-axis (>10°), and therefore we did not include 
these ISGRI data in our analysis of the persistent emission. 

For the long-term light curves we used OSA through the pro- 
duction of ISGRI images per single pointing in the 15-50keV 
range. This energy band was chosen to supplement the ISGRI 



^ |http: //integ ral ■ esac. esa. int/isda/| 

^ |http : //www . isdc ■ unige . ch/integral/download/osa.doc| 



light curves with the publicly available BAT long-term light 
curves (see Sect. I2.1.7l i. We force the flux extraction of each of 
the catalog sources, regardless of the detection significance of 
the source (see, e.g., Kuulkers et al. 2007). The correction for 
off-axis response is generally good up to ^10° from the centre 
of the field of view. We therefore selected only those pointings 
where 4U 0614+091 was less than 10° off axis. This led to a total 
of 1 15 pointings with different exposure times spanning the time 
interval from 2003 February 18 to 2006 April 18. 

2.1.7. 5w//f BAT 

The Burst Alert Telescope (BAT; Barthelmy et al. 2005) onboard 
Swift (launched November 2004; Gehrels et al. 2004) is a coded- 
aperture imager with a very wide field of view of about 2 steradi- 
ans, which operates in the 15-150keV band. The detector plane 
covers a net collecting area of 5200 cm^; the BAT angular reso- 
lution is 22' (FWHM). 

As soon as the BAT triggers on a y-ray burst. Swift auto- 
matically slews to the y-ray burst position and starts to observe 
the source with its more sensitive instruments. However, strong 
type I X-ray bursts may also trigger the BAT (but do not lead to 
automatic slews if they come from a known source). Assuming 
a black-body source with kT^2.5 keV, the 5cr detection limit for 
the BAT in 5 s is approximately 0.7 Crab in the 15-25 keV band. 

Two bursts from 4U 0614+091 triggered the BAT into a y- 
ray burst data collecting mode but did not result in automated 
slews. The triggers occurred on UT 2006 October 21 (nr. 234849 
on MJD 54029) and 2007 March 31 (nr. 273106 on MJD 54189). 
For these triggers, about 45 s of BAT event data were produced. 
A first account of these bursts has been given by Strohmayer et 
al. (2008). We ran the standard complete y-ray burst processing 
script on the BAT burst products of the two triggers to obtain a 
first impression of the burst behaviour, and to derive an updated 
position of the origin of the burst (see AppendixlAll. We then ap- 
plied the latest calibration available (as of August 2007) to these 
data (energy calibration, detector quality map, mask weighting), 
following the analysis threads provided by the Swift Science 
Centefl From these calibrated data we produced light curves, 
images and spectra. The first burst is strong enough to support a 
detailed time-resolved spectral analysis in the 15-30keV band, 
starting with a time resolution of 1 s. For the second, weaker 
burst we started with a time resolution of 2 s. We doubled the 
spectral time bin size, whenever the net-burst rates following the 
peak decreased by a factor of V2. 

The BAT continually monitors the sky; more than about 70% 
of the sky is observed on a daily basis. Results from this survey 
mode are publicly available in the form of light curves covering 
the 15-50keV energy band on two time scales: a single Swift 
pointing (^20 min) and the weighted average for each day0 In 
the latter a 6 mCrab source typically can be detected at the 3cr 
level (Krimm et al. 2006). We used the daily average light curve 
for 4U 0614+091 to study its long-term hard X-ray behaviour. 

2.2. Further data analysis issues 
2.2.1 . X-ray burst properties 

The band passes we used for the investigation of the burst light 
curves and extraction of some of the burst properties (see be- 



jhttp: //swift . gs£c .nasa. gov/docs/swift/analysis/ 
' ^http: //swi ft . gsfc. nasa. gov/docs/swift/results/transients/| 
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low) are 6-15keV, 2-28 keV, 1.5-12keV, 2-60 keV, 7^0 keV, 
and 15-30keV for the EURECAPNATCH, WFC, ASM, PCA, 
FREGATE, and both the ISGRI and BAT, respectively. 

For each of the ASM bursts, we derived the start time by 
taking the time when the count rate rose more than 5cr above the 
pre-burst level, while, for each FREGATE burst, the start time 
is the actual trigger time. For the other bursts we determined the 
start times by eye, because most of these bursts rise more slowly 
at higher energies. We define the burst rise time, t^^e, as the dif- 
ference between the time of the start and the time of the peak 
of the burst; the derived values were verified by visual inspec- 
tion of the light curves. We estimate the eiTor on the rise time to 
be about half a second (except during the spike of the burst ob- 
served on MJD51944, see Sect. 13.1. 2l i. The decay times, fdecay, 
are derived from fits of an exponential function plus a constant 
to the light curves. We fit from the maximum of the burst up to 
about 100 s after burst onset, if possible. For the longer bursts a 
longer time base was used when possible. When no post-burst 
information is available, we first fit the pre-burst rate with a con- 
stant and fix that when fitting the decay portion of the burst light 
curve. The range of values for which remains within 1 of the 
minimum value (i.e., Ax^ = l) is used to set the uncertainties of 
the decay times (assuming the fit parameters are uncoiTelated). 

Further burst properties include the bolometric fluence (Eb) 
and bolometric black-body peak flux (Fpeak), where the under- 
lying assumption is that the burst emission is on top of un- 
changed persistent emission (see, e.g., Kuulkers et al. 2002b), 
as well as information on the out-of-burst, i.e., persistent, emis- 
sion (see below). Estimates of E\, were obtained by integrating 
the fitted black-body models of the time-resolved spectra over 
photon energy and time (i.e., we do not take into account any 
residual burst emission outside the time range in which we did 
the time-resolved spectral fits). Since bursts decay more rapidly 
at higher energies, i.e., fdecay depends on the bandpass used, 
we also determine the characteristic decay time r=£'b/^peak 
which is not dependent on bandpass. Uncertainties in Eh and 
T were estimated by propagating the eiTors derived from the in- 
dividual time-resolved burst spectral fits and by assuming that 
these errors are symmetric. For the WATCH bursts and the 3 
weak bursts observed by FREGATE, time-resolved X-ray spec- 
tral analysis is not possible. Instead, we derive Fpeak from the 
peak count rates by scaling from the count rates and fluxes of 
the Crab. The Crab flux is taken to be 1.03x10"^ erg cm"^ s"' in 
the 6-15 keV WATCH band and 1.89x10"^ erg cm ^s ' in the 
7-40 keV FREGATE band. We then converted these values to 
estimated bolometric values assuming the emission comes from 
a 3 keV black body (which is fine for a large part of a burst seen 
at hard energies, see, e.g.. Fig. 1). For the 3 weak FREGATE 
bursts we derived Eb from the average burst-integrated spectra; 
for the WATCH bursts we derived E\, from the burst integrated 
number of counts and by performing a coiTection similar to that 
described above for Fpeak- We define the duration of the burst 
(fdur) as frise+2T, sincc not all bursts are fully covered. Inspection 
of the light curves show that these values are consistent with the 
total durations determined by eye. 

Burst oscillations at 414.7 Hz were present during a 5 s time 
interval in the brightest BAT burst cooling tail at the 4cr signifi- 
cance level (Strohmayer et al. 2008). We verified the presence of 
burst oscillations during this burst and searched for burst oscil- 
lations during the FREGATE and the 2005 ISGRI bursts where 
data were available at a high time resolution. We selected the 13- 
20keV band for BAT, 15-30keV for ISGRI and 7-40 keV for 
FREGATE, mainly to have the highest possible signal-to-noise 
ratio to search for a timing signature. We applied the Zj-statistic 



(standard Rayleigh statistic, see Buccheri et al. 1983) to the pho- 
ton event distributions for trial frequencies in a narrow win- 
dow centred on the expected frequency of 414.7 Hz (i.e., 413- 
416 Hz). The dynamical power spectra were computed using the 
same method as in Strohmayer et al. (2008). We employed in- 
tervals of 4 s or 8 s, stepping through the BAT, FREGATE and 
ISGRI bursts with intervals of 0.25 s; for the calculation of up- 
per limits we used intervals of 10 s. 

We searched data from both the PCA and HEXTE for os- 
cillations during the prompt burst seen on MJD51944 using 
fast Fourier transforms (FFTs). We used the single bit mode 
data available (full 2-60 keV PCA energy range taken at 122/is 
time resolution). FFTs of data from 4 s-long intervals within the 
burst were performed keeping the original time resolution. For 
HEXTE we computed FFTs of 1 s-long sets of 15^0keV band 
data with 122 fis time resolution. For the oscillation search in the 
long-lasting faint tails we employed the PCA data only. 

2.2.2. X-ray spectral fits 

In all our spectral fits (using XSPEC 11.3.0) we fix the low- 
energy absorption equivalent to a hydrogen column density of 
A^H=3xl02' cm-2 (e.g., Piraino et al. 1999; Mendez et al. 2002) 
for a composition of the absorbing matter as given by Wilms 
et al. (2000), employing absorption cross sections as provided 
by Verner et al. (1996). We derived errors in the fit parame- 
ters by using Ax^ - 1 (assuming the fit parameters are uncor- 
rected). We use a distance of 3 kpc (see Sect. 13.1.2b when re- 
porting black-body radii or emission areas. X-ray luminosities 
and absolute fluences. 

To describe the persistent emission we used ASM spectra 
around the ASM, FREGATE, 2003 ISGRI and BAT bursts, ASM 
plus ISGRI spectra around the 2005 ISGRI burst and WFC spec- 
tra around the WFC burst. We also used PCA-hHEXTE specti-a, 
whenever available within about 1 day of a burst. All spectra, 
except the PCA-hHEXTE spectra, could be well modelled by an 
absorbed power-law (power-law parameters: index, ypi and nor- 
malization at 1 keV, A^pi). For the PCA-hHEXTE spectral fits a 
soft component had to be included (we used a black body), as 
well as an FeK line in the form of a Gaussian (line parameters: 
line energy, Epe (keV), and photon flux, A^pe (photons cm"^ s"'); 
line width arbitrarily fixed at 0.1 keV). The time ranges over 
which the persistent ASM or WFC emission was integrated are 
given in Table[T] For the ASM spectra we generally used all data 
from about 1 .5 days before the burst up to the last dwell before 
the start of the burst. This leads to reasonably good spectra ex- 
cept for a few cases. 

To get an estimate of the bolometric (unabsorbed) persis- 
tent luminosity, we extrapolate the spectral fit results to the 
0. 1-200 keV band (arbitrarily chosen). Note that these values 
should be taken with some caution: power-law emission can con- 
tribute significantly to the total spectrum, at the lowest energies 
for high power-spectral indices (about 80% and 8% of the 0.1- 
200 keV emission is at 0.1-2keV and 10-200 keV, respectively, 
for ypi=2.5), and at the highest energies for low power-spectral 
indices (about 8% and 80%, respectively, for 7pi=1.5). 

XSPEC cannot provide errors on the integrated flux in an 
energy band. We therefore randomized the spectral parameters 
using the fit values and the derived Icr-eiTors; this was done 
10000 times and we recorded the resulting integrated 2-lOkeV 
and 0. 1-200 keV fluxes. The flux distributions are significantly 
skewed towards larger fluxes. We therefore fitted the flux dis- 
tributions below and above the peak of the distribution with 
Gaussians with different widths. 
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Fig. 1. Time profiles of raw full-bandpass photon rates (top panels), bolometric black-body flux (2nd panels), black-body tempera- 
ture (3rd panels), black-body radius for 3kpc distance (4th panels) and goodness of fit (bottom panels). The unit of intensity (top 
panels) is cts s"' for the 6 ASM bursts, cts s"' bin"' for the ISGRI burst, kcts s"' bin"' for the 2 BAT bursts, kcts s"' for the 10 
FREGATE bursts and cts s"' for the WFC burst. The bursts are denoted by instrument (FR=FREGATE) and MJD (see TablefB- For 
the ASM bursts which occurred on MJD 50200 and MJD 51 164 we show the Hght curves from SSCl and SSCO data, respectively 
(see Table©. 



Table 2. Burst and persistent emission properties for all bursts from 4U 06 14+09 1 . 
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"References: [1] Swank et al. (1978), [2] Mason et al. (1976), [3] Lewin (1976), [4] S. Kitamoto (priv. comm.), [5] Brandt et al. (1992), [6] http://space.mit.edu/HETE/Bursts/suinmaries.html 
* Superburst; " Also seen by //£7'£-2/WXM (H4087) [6]; ''Only seen by //£7'£-2/WXM [6]; " 2-20 keV; ^ 2-6 keV; « 2.5-7.5 keV; * 1.3-5 keV; ' 1-20 keV; ^ 6-150keV; * 6-12 keV; ' 20-100 keV. 
Note 1. The first 3 columns give the day of the burst, its start time (trigger numbers in brackets) and time since the last observed burst (A^)- For the burst parameters (fri.sc, fdecay, ^peak, E^, and fdm) see Sect. I2.2.T1 
X'^^/dof gives the goodness of the fit of the decay to an exponential and the degrees of freedom (dof). The rest of the columns refer to the average ASM persistent emission (Sect |2.2.'2l : time span over which it was averaged. 



goodness of fit of the emission to a power law and dof (X'^^ parameters y^i and A'pi. Fx,i, Fx and Fx.bol refer to literature" values of the absorbed persistent flux in various energy bands before 1994'' *, 

unabsorbed persistent fluxes at 2-lOkeV and 0. 1-200 keV after 1994, respectively. References are given between brackets, see footnote a. 
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The plus and minus errors in the integrated flux derived from the 
spectral fits were then taken as the Icr widths. 

2.3. Data and properties of previously reported bursts 

We complete the information on the properties of the bursts seen 
after 1992 with that of the bursts reported by Lewin (1976), 
Swank et al. (1978) and Brandt et al. (1992). We digitized the 2- 
60keV light curve of the 0S0-8/GCXSE burst from figure 2 in 
Swank et al. (1978) and also obtained further information on the 
event from that paper Data on the G/MA^Ar/WATCH burst (6- 
20keV; Brandt et al. 1992) is available in a form similar to that 
of the EURECA/WATCH bursts (see Sect. 12X1] ). We analysed 
this GRANAT/WAlTCH burst in a way similar to our analysis of 
the EURECAfWA^TCH bursts. Further information on the burst 
seen with the Horizontal Tube collimator (see Lewin et al. 1976) 
onboard SAS-3 was taken from Lewin (1976). Information about 
the persistent emission around the times of all these bursts, in- 
cluding those seen by EURECA/WAlTCH, were extracted from 
other literature (see Table[T]|. 

3. Results 

We report on 33 bursts from 4U 06 14+091, including the bursts 
previously reported and the superburst. They are listed in 
Table [1] Figs.[T](top panels) and|2]present the raw light curves 
(for a perspective on the short time-scale behaviour of the bursts) 
of, respectively, 20 unpublished bursts, and the 3 previously 
reported EURECA/WATCH and the 2003 ISGRl bursts. Time- 
resolved spectroscopy for the 20 bursts is also presented sepa- 
rately in Fig. [1] An overview of the normal bursts is given in 
Sect.[3ll] The burst that occurred on MJD 51944 and that shows 
evidence of radius expansion as well as a long-lasting faint tail is 
described in detail in Sect. 13.1.21 Our analysis of the superburst 
is presented in Sect. 13.21 

3. 1 . Overview of all normal bursts 

Apart from the fact that 4U 0614+091 is not bursting as sel- 
domly as previously thought. Table [T] shows that the recurrence 
time between bursts can be on the order of a week, when burst 
active (see also Sect. 13.5.1b . The bursts generally last between 
about 10 s and lOmin (ignoring long-lasting faint tails, see be- 
low). Based on the FREGATE bursts, we designate bursts as 
short bursts if the burst duration is less than about 100 s, and 
intermediate-duration bursts if they last between about 5 and 
lOmin. The intermediate-duration bursts thus refer to the two 
bursts that were seen by FREGAT E in 2 001 and 2002. The 
short bursts will be discussed in Sect. l3.1.Tl and the intermediate- 
duration bursts in Sect. 13.1.21 

Based on the facts that for all bursts the emission can be 
well described by black-body emission with canonical type 1 X- 
ray burst parameter values, that their profiles generally show an 
initial fast rise (at low energies) with a final slower decay which 
is exponential-like, and that the effective black-body temperature 
decreases during the decay for most of the bursts, we infer that 
the events are indeed type 1 X-ray bursts. 

3.1.1. Short bursts 

For the short bursts the rise times to maximum of the bursts as 
recorded by ASM, WFC and FREGATE are relatively fast, i.e., 
about 1-10 s. For die bursts seen by WATCH, ISGRI and BAT 



WATCH 49019 WATCH 49035 WATCH 49048 ISCRI 52867 




30 50 90 3D 60 90 30 60 90 30 60 90 



Time (sec) 

Fig. 2. From left to right: The 3 bursts observed with 
EURECA/WATCH in the 6-15 keV (top) and 15-150keV (bot- 
tom) bands, and the burst observed with ISGRJ in 2003 in 
the 15-30keV (top) and 30-200keV (bottom) bands. Zero 
time corresponds to the start of each burst (see Table [1]). For 
EURECA/WAlTCH detector count rates (in ctsbin ') are shown, 
i.e., no correction has been made for background and vignetting. 
The dashed line indicates the background level and the dotted 
line marks the 3cr excess level. For ISGRl the off-axis and back- 
ground corrected count rates (in cts s ') are shown. The bursts 
are denoted by instrument and MJD (see Table[T]). 



the rise times tend to be a bit longer, i.e., about 5-20 s; this can 
be attributed to these instruments being sensitive primarily at 
higher energies than the former instruments. Also, many of the 
bursts from these hard X-ray instruments show a two-stage rise, 
i.e., first a fast rise within seconds, then a slower rise to the peak 
intensity. The slow rise may be attributed to the slow increase 
in kT\,\, (see, e.g., Fig.lTJ, and therefore to the peak of the X-ray 
spectrum shifting towards the centre of the instrument bandpass. 

The short bursts have e-folding decay times between 2 and 
20 s; the shortest decay times are seen at the highest energies 
(ISGRJ and BAT), as expected if the burst spectrum softens dur- 
ing the decay. The characteristic decay time, r, is in the range 
6-42 s. 

The WFC burst (on MJD 51073) and two FREGATE bursts 
(on MJD 52961 and MJD 53740) are intrinsically weak; they 
reach peak intensities of less than a few Crab. The WFC 
event showed peak temperatures of about 1.5 keV, markedly 
below that seen during maximum of most of the other bursts 
(~3 ke V). These weak events may be similar to the burst reported 
from OSO-8, which reached a peak of ^0.36 Crab and showed 
/trbb=0.8±0.1keV during the first 20 s (Swank et al. 1978; see 
also Table[T]i. The event seen by SAS-3 (Lewin 1976) was rather 
weak too (see Table[TJ, but only part of a burst was seen, so one 
cannot infer its true peak intensity. 

The fluxes of the strongest bursts (including the 
intermediate-duration bursts) reach peak values of about 
15 Crab. Both this and the properties of the time-resolved spec- 
tral fits (see Fig. [T] and Table [T]) suggest they reach a limiting 
bolometric flux between about 21 and 32x10"^ erg cm"^ s"'. 
When taking into account only those bursts which reached a 
peak flux of more than 15x10"^ ergcm"^ s"', the average peak 
value is about 27x10"^ erg cm"^ s"' (fitting the observed values 
with a constant results, however, in a very high value ofx^^^, i-e.. 
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Fig. 3. Light curves (top row) and time-resolved spectral fit results (lower four rows) for selected bursts shown with higher time 
resolution than in Fig. [1] See the caption to Fig. [T]for brief explanations of the spectral parameters. The date (MJD) of each burst 
and the instrumental origin (FR=FREGATE) of the data are indicated. Left: FREGATE light curve and time-resolved spectroscopy 
for the burst on MJD 52997 with a time resolution of 0.16 s. The unit of intensity is kcts s '. Middle left: (top) PCA full band-pass 
Hght curve of the MJD 51944 burst with a time resolution of 0.001 s (from f =0-0.05 s), 0.01 s (from f=0.05-2.25 s), 0.125 s (from 
f=2.25-30s), and 16 s (from r=4000-8000s). The unit of intensity is lO^cs"' per 4 PCUs. The PCA (PCU3) rates after f=4000s 
have been multiplied by four. The ASM data are indicated with open squares; they have been scaled to the PCA rates using the Crab 
as a reference. Time f=0 s corresponds to UT 2001 Feb 4 21:52:42.8. (rest) Time-resolved spectral fit results using the data from the 
PCA (filled circles), ASM (open squares) and FREGATE (open triangles). Middle right: PCA full band-pass light curve (top) and 
time resolved spectroscopy of the long-lasting faint tail of a possible burst observed on MJD 51789. For this burst we do not have 
an exact start time, so f=0 s coiTesponds to the start of the data sequence (UT 2000 Sep 2 at 01 : 19:47). Intensity is in units of c s"' 
for PCUs 2 and 3 combined. Right: ASM full-bandpass Hght curve of the superburst on MJD 53441 from 4U 0614-1-091 (top) and 
time resolved spectroscopy. Time is given in hours since UT 2005 March 12, Oh, and the intensity is given in cts s"' . 



23, for 14 do£3)- The highest bolometric peak flux measured is 
31.5+0.5x10"^ erg cm^^s"'. 

The observed fluences, E\,, of the intrinsically weak bursts 
are between 7 and 13x10"^ erg cm"^. The values for the other 
short bursts range from about 16 to about 124x10"^ erg cm"^. 
For a source distance of 3 kpc we derive absolute fluences from 
the above quoted values between 0.8 and 13x10^'^ erg. 



The goodness of fit is expressed by the reduced chi-squared, or;^'^^^. 
It is the sum of the weighted square deviations between the data and the 
model, divided by the number of degrees of freedom (dof). 



We also inspected the burst light curves at higher time resolu- 
tions. One of the FREGATE bursts (MJD 52997) showed an in- 
teresting spike at the start of the burst. The high-time resolution 
light curve is shown in Fig.[3](top left). The spike lasts for about 
0.3 s; after that the flux slowly rises to maximum for about 4 s. 
Double peaked light curves are typically seen in hard X-ray light 
curves when photospheric radius expansion is important (see, 
e.g., Lewin et al. 1984, Tawara et al. 1984, Cocchi et al. 2000). 
We performed time-resolved X-ray spectral fits at higher time 
resolution, to follow in more detail the evolution of the burst, es- 
pecially the spike (Fig.|3] left, bottom panels). During the spike. 
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Fig. 4. PC A 2-60 ke V light curve centred around the burst onset 
on MJD51944 at a 10 msec time resolution. The intensity is as 
measured by 4 PCUs combined. Time t-Q s corresponds to UT 
2001 Feb 4 21:52:42.8. 

the spectrum is quite hard compared to that generally observed 
near maximum of type 1 X-ray bursts (^rbb=4.3;^y jkeV versus 
^3keV), and R\,b is quite small. During the first few seconds 
Rbb rises, after which it drops slightly and then levels off. kT^b 
quickly drops after the spike, but then it increases up to ^3 keV 
at the peak of the burst. This may be a hint of a short radius- 
expansion phase, but the statistics are somewhat poor. Note that 
Fboi does not stay constant but slowly rises during that phase, 
casting some doubt on the radius-expansion interpretation. 

3.1.2. Intermediate-duration bursts and evidence of 
pliotosplieric radius expansion 

The PCA light curve around the time of the start of burst 
MJD 51944 is shown in Fig.|4] One may see a strong, very short, 
spike followed by a large drop in the flux, and then the onset 
of the main burst. In the top part of Fig. |3] (middle left), the 
first part of the burst is plotted again, but the spike is shown 
at higher time resolution. During the spike the intensity reaches 
maximum within 5 ms. Then, within 0.04 s, the intensity drops 
to a level consistent with the background and certainly below 
that seen before the spike. About 1.1s later emission reappears 
again and continues as the main bursO- A local maximum is 
then reached within 2 s. In the spike, higher count rates are ob- 
served than during the first 30 s of the main burst. 

The spike, drop below pre-burst levels and main burst in- 
dicate strong radius-expansion (see, e.g., Molkov et al. 2000, 
Strohmayer & Brown 2002; see also in 't Zand & Weinberg, in 
preparation). The results of the PCA time-resolved spectral fits 
are shown in the bottom panels of Fig. |3] (middle left), together 
with spectral fits to the average of the ASM measurements about 
0.3 hrs after the start of the burst. We also overlaid the spectral fit 
results from FREGATE (see Fig. 1, most left of bottom panels). 
Indeed, the inferred emitting area clearly increases during the 
spike and decreases during the very first part of the main burst. 
The expansion is more than 3 orders of magnitude. The peak of 
the emission (as parametrized by kT^b) therefore shifts outside 

' ' FREGATE triggered on the main burst. 



the PCA observable energy band, causing the dramatic drop in 
intensity. 

After the fast contraction, it takes about 100 s for the burst 
to reach a maximum inferred temperature. During that time the 
emission area slowly decreases further (the end of this phase is 
usually referred to as 'touch-down'). During the maximum of 
the main burst, between about 5 and 100 s after the start, the flux 
varies only slightly between about 18 and 26x10"^ erg s"' cm"^. 
This is consistent with the peak flux values reached during the 
other strong bursts (see Sect. 13.1.1b . If we assume that during 
that time the flux is close to the Eddington value, the distance 
to the source may be inferred. By using the highest measured 
peak flux (Sect. 13.1. Il l and the empirically derived Eddington 
luminosity of 3.79x10^*^ erg s ' appropriate for a pure He accre- 
tor (Kuulkers et al. 2003), we infer a distance to the source of 
3.2 kpc. Since the estimated uncertainty in the empirically de- 
rived Eddington luminosity is about 15% (see Kuulkers et al. 
2003), we use for simplicity a distance value of 3 kpc in all our 
calculations. We note that during the strong expansion and con- 
traction phase the apparent flux increases by more than an order 
of magnitude with respect to that described above. This would 
mean that super-Eddington fluxes are reached; however, since 
most of the emission is outside the PCA sensitive energy band 
and the spectra in this phase are highly non-Planckian (see be- 
low), these flux values must be regarded with caution. 

Because of the very high count rates, the spectra constructed 
from the PCA data have small statistical errors, and this of- 
ten means that simple spectral models are not sufficient to fit 
the spectra within the uncertainties. We find;^^^^^^ values as high 
as 20, especially during the first part of the contraction phase 
(see bottom panel of Fig. |3] middle left). Deviations from a 
pure black body are seen frequently during radius-expansion 
bursts (see, e.g., Kuulkers et al. 2002b, 2003, and references 
therein). Since the temperature changes are fast during the ex- 
pansion/contraction phase, a single black body may not be ade- 
quate in the time bins used. However, 10-20 s after the start of 
the bursts the temperature stays more or less constant, while still 
the spectra fits are not ideal. The deviations are thus intrinsic to 
the source. 

Both intermediate-duration bursts show very long rise times 
during the main burst part, i.e., ^70 s and ^ 140 s, on MJD 5 1944 
and MJD 52322, respectively. They also have long exponential 
and characteristic decay times, i.e., fdecay-40 s and ^90 s, respec- 
tively, and T^l 19 s and ^216 s, respectively.. 

The burst on MJD 5 1944 showed a long and faint tail, last- 
ing for hours: after the main burst, during the next pass of PCA 
observations (between about 1.2 and 2.2 hrs after the start of 
the burst), the intensity was still slowly decaying but had not 
yet reached pre-burst levels. Serendipitously, the ASM scanned 
over 4U 0614+091 about 0.3 hrs (i.e., in between the two above 
described PCA passes) and 1.9 hrs after the start of the burst. 
The trend seen in the ASM data points is consistent with that 
seen in the second PCA pass (see Fig. [3] top middle left). The e- 
folding decay as measured from the bolometric black-body flux, 
Fboh light curves, Texp,bb, is 1 143;^4g2'^ s, i.e., a factor of about 35 
longer than that measured during the first part of the decline from 
maximum of the main burst (Texp,bb=33^^ s). During the first part 
of the decline kT^b drops quickly (see Fig.[T]), but after the flux 
has dropped by more than an order of magnitude, kTbb does not 
decrease anymore, but stays more or less constant near 1 keV; the 
decrease in flux is then due to an inferred decrease in emitting 
area. About 5.2 hrs after the start of the burst the flux still had not 
yet reached pre-burst levels; the steady flux decay suggests that 
it is still due to remaining burst emission. However, at that low 
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Fig. 5. HEXTE hard X-ray (15-60keV, cluster A) light curve of 
the burst on MJD51944of 4U0614-H091 (filled circles), as well 
as the superburst light curves of 4U 1820-303 (filled triangles) 
and 4U 1636-536 (open squares). The time resolution is 16 s 
for 4U 0614-1-091 and 32 s for 4U 1820-303 and 4U 1636-536. 
Time = Os sec corresponds to UT 2001-02-04 21:52:43, 1999- 
09-09 01:46:54 and 2001-02-22 16:52:11, for 4U0614-H091, 
4U 1820-303 and 4U 1636-536, respectively. 



flux level we cannot be certain whether this is indeed remaining 
burst emission or due to a slight change in the persistent flux. 

We found one other clear example of a long and faint tail 
in the PCA data, on MJD 51789, between UT 01:19:34 and 
02:10:35 (see Fig. [3j top, middle right). In the bottom pan- 
els of Fig. [3] (middle right) we show the results of the time- 
resolved fits. The behaviour is similar to that seen in the long- 
lasting faint tail of the burst on MJD 5 1944, albeit at a sHghtly 
higher flux level: kT^b stays constant near ^1 keV, while the in- 
feiTed emitting area decreases by a factor of about 2. The ex- 
ponential decay time as measured from the bolometric black- 
body flux light curves is Texp,bb=868+32s. No RXTE observa- 
tions of 4U 0614+091 were obtained for about 5 hours prior to 
the start of the observation at UT 01:19:47 and none for about 
2.5 hours following the end of this observation. Therefore, an 
intermediate-duration burst like that seen on MJD 5 1944 may 
have preceded the tail here. 

Although the long-lasting faint tails last as long as vari- 
ous superbursts, their properties are different. Superbursts still 
show evidence of cooling, hours after they started. We show 
this in Fig. |5]by comparing the hard X-ray (15-60keV) emis- 
sion during the intermediate-duration burst with that seen during 
the two superbursts observed from 4U 1820-303 (Strohmayer 
& Brown 2002) and 4U 1636-536 (Strohmayer & Markwardt 
2002). Clearly, hard X-ray persists during most of the superburst, 
while it only lasts for minutes during the intermediate-duration 
burst. 

For the two intermediate-duration bursts the observed flu- 
ences are approximately 320x10"^ erg cm"^ (MJD 5 1944) and 
610xl0"^ergcm"^ (MJD 52322), corresponding to absolute flu- 
ences of 34x10"'^ erg s"' and 67x10-'^ erg s"', respectively, at 
3 kpc. We have ignored here any long-lasting faint tails; we 



expect that corrections for contributions from the tails would 
change the above values by only a few per cent (see Table[T]i- 

3.2. Superburst 

Inspection of the ASM light curve of 4U 0614+091 revealed 
a long flare on March 12, 2005 (MJD 53441; see top right of 
Fig. |3] and Fig. |6] see also Kuulkers 2005). Between about UT 
10:38 and 16:53 the ASM flux increased by a factor of 9.5 up 
to 0.3 Crab. About 1.5 hours later the flux had dropped to 0.17 
Crab; ^9.5 hours later it was still a factor of about 2 above the 
pre-flare flux level. The exponential decay time, tjecay, of the 
flare is ^2.1 hr. We calculated the hardness ratio using various 
energy bands. We found the most significant results when we 
used the ratio of count rates in the 3-12keV to 1.5-3 keV bands, 
as shown in Fig.|6] The hardness was higher during the flare than 
in the pre- or post-flare measurements; it softened during the de- 
cay. We investigated the raw dwell light curves during the flare, 
and found no evidence of strong variability, anomalous features 
such as normal bursts, or for instrumental artifacts. We, there- 
fore, confidently attribute the high fluxes to the source. 

We grouped the flare into three time bins, and performed a 
time-resolved spectral analysis of the average ASM spectra in 
these time bins. The spectra of the first two time bins are in- 
consistent with an absorbed power law (x^^Jdof of 6.1/1 and 
16.2/1, respectively; for the third time bin ;i'^^ydof =0.7/1). The 
net-flare spectra can best be modelled by absorbed black-body 
spectra. The results of this time resolved spectral analysis are 
displayed in Fig. [3] (right). The peak flux is rather low, about 
0.7x10"** ergs"' cm"^, and the maximum value for kT^b black- 
body temperature only about 1 .5 keV. However, we may not have 
seen the true peak intensity or black-body temperature due to 
the sparse coverage by the ASM. The inferred emission areas 
are rather low (2^ km). They are similar to those inferred for 
the long-lasting faint tails of the bursts observed by the PCA 
(Sect. 13.1.2b . However, kT^b is somewhat higher, and the time 
scales involved about a factor of 7 longer, than those seen in 
these tails. Note that the inferred radii are not inconsistent with 
those seen during some of the superbursts from other LMXB 
burst sources (e.g., Kuulkers et al. 2002a). Since our observed 
flare resembles other superbursts, we suggest the flare to be a 
superburst. 

3.3. Search for burst oscillations 

We confirm the presence of burst oscillations during the bright- 
est BAT burst, but do not detect any significant signal in any of 
the other bursts. This is no surprise: compared to the bright burst 
BAT light curve, the 2005 ISGRI burst light curve has a fac- 
tor 6 lower signal-to-noise ratio, and the FREGATE light curves 
have a factor 1 .6 to 2 lower signal-to-noise ratio. Therefore, a 4cr 
detection by BAT would coiTespond to non-detection by both 
ISGRI and FREGATE, assuming that the fractional amplitude 
was equal to that measured in the BAT data. We derive 3cr upper 
limits on the fractional rms (e.g., van der Klis 1995) in the 413- 
416 Hz range of about 4-6% at the peaks of the bursts observed 
by FREGATE, while they are about 9-12% in the burst tails. 

The 3cr upper limits on the fractional rms amplitude of burst 
oscillations during the prompt burst on MJD 5 1944 are about 
0.3% and 2% in the 100-1000 Hz range as measured by the PCA 
and HEXTE, respectively. Thus, if there were oscillations during 
the prompt burst on MJD 5 1944 they were much weaker than 
those observed during the brightest BAT burst on MJD 54029. 
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Table 3. Results of spectral fits of PCA+HEXTE data taken outside of but within one day of bursts, i.e., representing persistent 
emission only.^ 



MJD 


Af 






^bb 


^bb.bol (10 


rpi 


A'pi 




A^Fe (10-^ 


Fx 


^X.bol 




(days) 




(keV) 


(km) 


erg s"' cm"-) 






(keV) 


ph s"' cm"^) 


(lO^'^erg 


s cm -) 


51789 


<0.27 


1.81/108 0.93±0.01 


1.45±0.01 


1.92+0.03 


1.841+0.003 


0.206+0.001 


6.9+0.1 


0.52+0.05 


8.37+0.04 


35.9+0.2 


51944 


0.09 


1.17/108 


1.02±0.02 


1.18±0.01 


1.75+0.04 


1.988+0.004 


0.289+0.001 


6.6+0.1 


0.58+0.08 


9.10+0.06 


37.7+0.3 


53040 


1.05 


1.14/108 


1.09±0.01 


1.12+0.01 


2.19+0.04 


2.162+0.004 


0.269+0.002 


6.8+0.2 


0.34+0.11 


7.36+0.06 


29.7+0.2 


54188 


1.02 


1.35/79" 


0.90±0.02 


1.51+0.02 


1.86+0.04 


1.957+0.003 


0.256+0.001 


6.7+0.1 


0.68+0.07 


8.44+0.05 


35.3+0.2 



" No HEXTE Cluster A spectrum. 



Note 2. Given are the day of the observation (MJD), the time to the closest burst (Af), and the spectral parameter values. The latter are the goodness 
of fit degrees of freedom {Y^,^Jdof), the black-body parameters (temperature kTf,f„ radius i?bb at 3kpc, bolometric black-body flux Fbb,boi)> the 
power-law parameters (power-law index ypi, normalization at 1 keV A'pi), the Gaussian line (presumably Fe K) parameters (line energy Ep^, line 
normalization A'pe; line width was fixed at 0.1 keV), and the unabsorbed fluxes in the 2-lOkeV (Fx) and 0.1-200 keV (Fx_boi) bands. 
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Fig. 6. ASM individual dwell intensity measurements (1.5- 
12keV) and average hardness ratios around the time of the su- 
perburst from 4U0614H-091. Hardness is defined as the ratio of 
the intensity in the 3-12 keV band to that in the 1 .5-3 keV band. 



No evidence of oscillations was found during the faint long- 
lasting tails either, with 3cr upper limits of about 1.2%. 

3.4. Persistent emission 

Moderate resolution X-ray spectra of the persistent emission 
from 4U 0614-1-091 above typically 1 keV have been satisfac- 
torily fit by a soft component plus a hard, power-law compo- 
nent, both subjected to interstellar absorption. Generally, for 
4U 06 14-1-091 photon indices for the power-law component are 
found to be between 2 and 3 and temperatures for the soft black- 
body component are between kT -0.5 and 1.5 keV. Total fluxes 
(1-20 keV) are typically between 1 and 5x10""^ erg s"' cm"^ 
(Barret & Grindlay 1995, Ford et al. 1996, 1997; see also 
Piraino et al. 1999, Fiocchi et al. 2008). Our four most accu- 
rate measurements for these components are provided in Table|2l 
Our other lower-quality measurements based on power-law fits 
only are provided in Table [1] Our derived values of the power- 
law index, F, vary between ^1.6 and ^2.5. The persistent 2- 
lOkeV flux is around 0.5-2x10"^ erg s"' cm"^, with an average 
of 0.66+0.06xl0"'erg s ' cm"^. Our spectral fits are consistent 
with that derived previously. 



The unabsorbed bolometric flux (as estimated from the 0.1- 
200 keV flux estimates in Table [1] see Sect. 12.2.2b is about 
a factor of 5 higher than the 2-lOkeV flux, i.e., on average 
3.4+0.2x10""^ ergs"' cm"^. The best estimates of the bolomet- 
ric correction factor from the 2-lOkeV flux are obtained from 
a joint analysis of the PCA and HEXTE spectra, because of the 
quality and broad-band coverage. They are between 4.0 and 4.3 
with an average of 4.16+0.05 (see Table|2]i, which is consistent 
with that derived from the ASM plus ISGRI spectra (4.4+1.3, 
see Table [1]). These factors are compatible with that derived by 
in 't Zand et al. (2007) from a broad sample of LMXBs, i.e., 
a factor 2.9+1.4 between the 2-lOkeV flux and the bolometric 
flux as estimated from the extrapolated 0. 1-100 keV flux, but 
are a bit higher than the value estimated by Migliari & Fender 
(2006) for an atoll source in the hard state, i.e., 2.5. Our bolo- 
metric fluxes derived from the ASM spectral fits are consistent 
with the fluxes derived from the broad-band PCA plus HEXTE 
spectral fits. Also, our fluxes are comparable to that measured 
using broad-band BeppoSAX dsitSL on 4U 0614+091 in the 0.1- 
200 keV band (Piraino et al. 1999; see also Fiocchi et al. 2008), 
as well as simultaneous Swift/XRT and PCA plus HEXTE data 
in the 0.6-100 keV band (MigHari et al. 2009). This lends confi- 
dence to the bolometric correction and the interpretation that the 
extrapolated 0. 1-200 keV flux represents the bolometric flux. 

If the source emits isotropically during and between bursts, 
and if the peak flux during the bursts Fpeak represents emission at 
the Eddington limit, then the parameter y (=/^pers/^peak, where 
both fluxes are bolometric) gives the source persistent luminos- 
ity in terms of the Eddington limit (see, e.g., Cornelisse et al. 
2002b). Using the values from Table [T] and using only those 
bursts where Fpeak>15xlO"^ergs"' cm"^ (see Sect. B.l.TT i. we 
find values for y between about 0.0032 and 0.016, with an aver- 
age of y=0.0 13 ±0.001. This is consistent with the estimates by, 
e.g.. Ford et al. (2000) and in 't Zand et al. (2007). 

3.5. Long-term X-ray beliaviour 

We investigated the long-term X-ray behaviour of 4U 0614+091 
to see whether it affected its bursting behaviour using results 
from the RXTE/ ASM (Sect. 13.5.1b and instruments on a number 
of other spacecrafts (Sect. l3.5!2l i. 



3.5.1. ASM long-term persistent flux and burst recurrence 
times 

The ASM data provide a uniquely long and homogeneous his- 
tory of the persistent flux. In Fig. |7] (left) we show the ASM 
(1.5-12keV) long-term light curves as daily averages. The cor- 
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Fig. 7. Left: Daily averaged ASM light curve (1.5-12keV). The times of the normal bursts (see TablelTJ are indicated with arrows 
at the top. The superburst occurrence is marked with an 'S'. Right: The yearly averaged ASM light curve (1.5-12keV; filled circles) 
and the measure of the variability (rms) in the yearly time bins (filled triangles). 

Table 4. Average burst recurrence time estimates for various instruments during three periods. 



1996-2000 flaring period (I) 2001-2002 calm period (II) 2003-2007 flaring period (III) 



Inst 


^exp 


n" 


Ab" 


Inst 


^exp 


n 


6tb 


Inst 


^exp 


n 






(d) 




(d) 




(d) 




(d) 




(d) 




(d) 


WFC 


24.9 


1 


11-83 


WFC 


2.3 





>5 


WFC 








ASM 


27.5 


2 


7-31 


ASM 


9.6 





>19 


ASM 


20.5 


4 


3-9 


PCA 


9.6 





>19 


PCA 


3.7 


1 


2-12 


PCA 


9.4 





>19 


FREGATE 








FREGATE 


^65 


2 


17-72 


FREGATE 


^129 


11 


9-16 


ISGRI 








ISGRI 








ISGRI 


25.3 


2 


7-28 


All 


^62 


3 


12-39 


All 


==81 


2' 


22-90 


All 


^184 


17 


9-14 



" n = number of bursts observed in the considered period. 

' (5?b = expected average burst recurrence time based on Poisson likelihood Icr limits, see text. 
^ Note that PCA and FREGATE observed one of the bursts simultaneously, see Sect. 13.1.21 



responding average flux and variability (rms) per year are shown 
in Fig. Q (right). From 1996 to mid 2000 the source was flar- 
ing on time scales of days to weeks from about 25-50 mCrab 
up to about 100-125 mCrab, with yearly averages above 30- 
35 mCrab. Then the source stayed at a relatively quiet level near 
30mCrab for about two and a half years (mid 2000 to 2003). 
Flaring behaviour similar to before mid 2000 was seen from 
2003 to mid 2005. From mid 2005 to mid 2007 4U 0614+091 
was still flaring, but less intensely, i.e., it varied from 25- 
50mCrab up to about 75-lOOmCrab and the yearly averages 
were again above 30-35 mCrab. The figure suggests that the 
degree of flaring and the long-term average (months to years) 
intensity seem to be related; a high degree of flaring is seen 
when the yearly average persistent flux is above about 30- 
35 mCrab and a low degree of flaring is seen otherwise. A flux 
of 30-35 mCrab (1.5-12keV) coiTesponds to about 1% of the 
Eddington limit (see Sect. 13.41 ). For the following discussion we 
roughly divide the ASM history of 4U 06 14+091 into three pe- 
riods: (I) 1996-2000 flaring period, (II) 2001-2002 calm period, 
and (III) 2003-2007 flaring period. 

In Fig. Q (left) we mark the times of the normal bursts, as 
well as the superburst. The timing of burst detections is related, 
in part, to the times of year when the Sun is not near the source 
and the latter is more easily observed. Moreover, bursts could 
have occurred during many of the data gaps. It is also possi- 
ble that with the ASM we missed the bursts which last signifi- 
cantly longer than the dwell duration (like the two intermediate- 
duration bursts). Also, weak bursts like that seen by the WFC 



may have been gone unnoticed with the other instruments; more- 
over, FREGATE, ISGRI and BAT operate in a higher bandpass 
than the ASM and WFC, and have to cope with higher back- 
grounds and so tend to trigger on brighter bursts. Therefore, we 
cannot rule out that we did not see all bursts which happened in 
our time frame of interest. 

Fig- 13 (left) suggests, however, that bursts appear more fre- 
quently during the 2003-2007 flaring period with respect to the 
previous periods. This is especially apparent if we take into ac- 
count that during the 2000-2002 period there were 3 surveying 
X-ray instruments active (WFC, ASM and FREGATE). To sup- 
port this suggestion, we determined the average burst recurrence 
times during the three above described periods for each instru- 
ment, as well as all together. We excluded results from the BAT 
(we only have information on triggered bursts, i.e., we did not 
perform a systematic search through the BAT archive) and JEM- 
X (very low source exposure time). We also did not take into 
account either the tail of a possible burst on MJD 51789 or the 
superburst. When considering all instruments, we assumed that 
the exposure intervals did not overlap, and that they have the 
same efficiency in detecting bursts. We also assumed that the ex- 
posure time of a particular instrument is randomly distributed in 
each of the three periods. If n is the number of bursts seen during 
a certain period, then the average burst recurrence time is defined 
as fexp divided by n, where fexp is the total exposure time on the 
source in question in that period (note that this does not necessar- 
ily mean that the bursts recur periodically). Since we are dealing 
with relatively low number of observed bursts, Poisson statis- 
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tics apply. For the expected range in average burst recurrence 
times ((5fb) we use the Poisson likelihood Icr lower and upper 
limits on the expected number of events when n events are ob- 
served (where n=0,l,2,3,- ■ ■)■ We verified our method by doing 
Monte Carlo simulations. The net source exposure times, num- 
ber of bursts observed, and resulting average burst recurrence 
times per period of interest are shown in Table |4] For the indi- 
vidual instruments, the ASM and FREGATE give the best con- 
straints on the average burst recurrence times, mainly because of 
the number of bursts observed. The estimates of the PCA aver- 
age burst recuiTence times differ significantly from those for the 
other instruments; we attribute this is to the fact that the data are 
not taken serendipitously, as well as to the relatively low expo- 
sure times. 

The average burst recurrence times in periods I and II are not 
significantly different; the average over the two periods com- 
bined is 29^[q days (taking into account all considered instru- 
ments). Period III shows a significant change in the average burst 
recurrence time with respect to the previous period: 11 days 
(taking into account all considered instruments); it is just con- 
sistent with period I. The lower limit on the average burst re- 
cuiTence time in period III is close to the shortest observed 
value of At in that period (7 days, see Table [TJ. During the 
EUREC A/WATCH period in 1993 the burst average recurrence 
time is 17^ days, consistent with values of At in that period 
(see Table [Ij. It is consistent also with the values of all three 
ASM periods I, II, and III. 

Interestingly, the only two bursts that were observed in 2001 
and 2002 are both intermediate-duration bursts. They thus oc- 
curred during a period when 4U0614H-091 was not only less 
burst active, but also when the persistent emission was rather 
calm and had the lowest averaged values. All other bursts (22, 
including the BAT bursts) are shorter and occurred during the 
two flaring periods. The burst duration, therefore, seems to be 
related to the level of the persistent emission, and possibly to the 
degree of flaring. 



3.5.2. Further long-term broad-band X-ray behaviour 

4U 0614-1-091 was discovered by UHURU with variable flux lev- 
els between about 15 and 70mCrab (2-6 keV; Giacconi et al. 
1972). Other early measurements of the flux, besides those from 
UHURU, may be found in, e.g.. Mason et al. (1976), Swank et 
al. (1978), Parsignault & Grindlay (1978), Markert et al. (1979), 
and Warwick et al. (1981). In that time it was always seen at a 
flux of about 10-lOOmCrab in roughly the ~2 to ~10keV band. 
During sporadic measurements by the G/«^fl/All Sky Monitor 
between February 1987 and November 1991 4U 0614+091 was 
never detected in the 1-20 keV band, with an upper limit on the 
source flux of about lOOmCrab (S. Kitamoto, private communi- 
cation). 

4U 0614+091 was too weak to be detected by 
EURECA/WATCH on a daily basis, but data from the en- 
tire observation period can be combined. By adding the 
skymaps (see Brandt 1994 for a description of the technique), 
4U 0614+091 was detected at an average persistent flux of 
25+lOmCrab in the 6-12 keV energy band (see also Fig. lA.ll 
right). 

We investigated the CGRO/BAlTSE occultation measure- 
ments (see, e.g., Harmon et al. 2004) in the period April 
1991 and November 2000; they show 10-day average 20- 
lOOkeV fluxes from non-detections (with typical 10-day 3cr up- 
per limits between about 50-70 mCrab) up to about 150mCrab. 
During the period of the EURECAfWALTCH observations when 
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Fig. 8. Zoom in on the long-term hard X-ray (15-50keV) light 
curve (top) and soft X-ray (1.5-12keV) light curve (bottom), 
of 4U 0614+091. Shown are the 1-day and 10-day averages for 
ISGRI (filled circles) and Swift/BAT (grey filled triangles), re- 
spectively (top), and 10-day averages for the ASM. Time in 
years is given above the plot. The occuiTence times of the normal 
bursts are indicated by aiTows at the top. The superburst occur- 
rence is marked with an 'S'. 



4U 0614+091 was in the field of view, the 20-100 keV flux was 
on average 57+6mCrab; close to the time of the bursts seen by 
EURECA/WATCn it was ai'ound lOOmCrab. 

The BAT and ISGRI 15-50keV long-term light curves are 
consistent with each other. The 15-50keV flux in the February 
2003 to September 2007 time frame ranges from a few mCrab to 
50 mCrab, but is most frequently around 20 mCrab. A zoom in 
on the 15-50keV light curve and the corresponding part of the 
ASM light curve is shown in Fig. |8] Three interesting features 
may be seen in this figure. First, there is a rapid increase to about 
50 mCrab for about a month, followed by a drop within weeks 
to about lOmCrab at MJD 53500. After about 2 months it goes 
back up again to the average level. Second, whenever the hard 
X-ray flux drops below roughly lOmCrab, the soft X-ray inten- 
sity increases, up to a factor of 2 or so. Such an anti-coiTelation 
between soft and hard X-rays on these long time scales was pre- 
viously discussed by Barret & Grindlay (1995) and Ford et al. 
(1996). On the other hand, hard X-ray flares are not accompa- 
nied by changes in the soft X-ray emission. Third, the super- 
burst occuiTed around the time when the 15-50keV flux was 
at its highest level (unfortunately no simultaneous BAT mea- 
surements exist at the time of the superburst). Just after the su- 
perburst the 10-day average 15-50keV flux is low, while the 
1.5-12keV flux is at high levels. The ASM data show that this 
high average flux is not only due to the superburst itself, but 
also to a somewhat increased persistent level just after the su- 
perburst (see Fig.|6]l: before the superburst the 1.5-12keV flux 
was 31 + 1 mCrab, whereas for several days after the superburst 
it was at an elevated level of 53+2mCrab. About 10 days after 
the superburst it was back at the pre-superburst level. A normal 
burst (MJD 53460) was seen only 19 days after the superburst 
(see Table [TJ; note that this is only slightly longer than the typ- 
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Soft Colour Intensity (mCrab) 

Fig. 9. PCA colour-colour diagram (Left) and hardness-intensity 
diagram (Right) of all public data of 4U 06 14+091. Plotted are 
all persistent data points averaged per observation, and nor- 
malised using the Crab (see Sect. 12.1.3b . Data taken close to a 
burst are marked with a black filled triangle (Af<0.3 days) or a 
black filled circle (Af~l day). For the definitions of the colours 
and intensity, see Sect. 12.1.31 



ical average burst recurrence time when 4U 06 14+091 is burst 
active (see Sect. 13.5b . 

The CD and HID (Fig. |9]l are qualitatively similar to those 
presented for the PCA data obtained in 1996-1998 by van 
Straaten et al. (2000). The colours often change on a daily time 
scale, as is typical for atoll sources (e.g., Hasinger & van der 
Klis 1989). For a wide range in hard colour, the soft colour does 
not change noticeably. Only near the lowest hard colours does 
the soft colour start to change; then there is an anti-correlation 
between hard and soft colour When the hard colours are low- 
est, however, the soft colour changes are large while the hard 
colour changes are small. This is in line with the broad-band 
behaviour on longer time scales of the fluxes in the 1.5-12keV 
and 15-50keV bands described in the previous paragraph. In the 
HID higher intensities correspond to lower hard colours. Again, 
at the highest intensities (above about 50 mCrab), however, the 
hard colour changes only very little. The part of the diagram 
where the soft colour and intensity increases at constant hard 
colour (value of about 0.5) is referred to as the soft state or 'ba- 
nana' branch (see van Straaten et al. 2000). The rest is referred 
to as the intermediate state or 'island' branch. 

For 4 bursts (i.e., those that occurred on MJD 51789, 
MJD 51944, 53040 and 54188) we have information about the 
persistent emission within about 1 day of the burst (see Table|2]i. 
In Fig.|9]we indicate the position the persistent source was in at 
these times. All 4 bursts occurred when the persistent emission 
was hardest and when the intensity was below ^40mCrab, i.e., 
when the source was in the intermediate or hard state ('island' 
or extreme 'island'). 



4. Discussion 

4.1. Observational summary 

We have expanded the number of observed bursts that can be at- 
tributed to 4U 0614+091 from a handful to 33. These include a 
superburst. Most of the newly-found bursts are localized within 
1'.5 of the optical position of 4U 06 14+091. This confirms the 
type I X-ray burster nature of this low-mass X-ray binary. We 
find a significant change in the average burst recurrence time in 
a decade of observations, i.e., from roughly a month before mid- 
2003 to about one to two weeks afterwards. When burst active, 
the bursts have peak luminosities of about 3x10^** erg s"', last 



for up to about a minute, and release energies of about lO-'''- 
10"^" erg. The larger set of bursts shows a wide variety of charac- 
teristics, with dynamic ranges of a factor of 40 in peak flux, 100 
in duration, and 230 in absolute fluence. 

We found three long-duration bursts. Two of them have ini- 
tial decay time scales of about 100 s and resemble intermediate- 
duration bursts seen from other sources (thought to be due to 
deep He flashes; see in 't Zand et al. 2005 and Cumming et al. 
2006). The third has an initial decay time of about 8000 s, com- 
parable to the decay times of superbursts (due to an even deeper 
C flash contained in the ocean of the NS envelope; see Cumming 
& Bildsten 2001 and Strohmayer & Brown 2002). It lasts a fac- 
tor of 10 longer than the longest known intermediate-duration 
burst (for example, see in 't Zand et al. 2007, Keek & in 't Zand 
2009). Unfortunately, the coverage of this burst is incomplete. 
This prevents a measurement of the true peak flux and a check 
for the occurrence of radius expansion. The average burst re- 
currence time when the two intermediate-duration bursts were 
seen, is notably longer than when the shorter bursts were seen 
(40;*^jg days during a two-year period versus 1 l^j days for other 
bursts in the subsequent 4.7-year period). 

One of the intermediate-duration bursts shows clear signs 
of strong radius expansion within the first few seconds of the 
event. The flux seen at the end of this episode of radius expan- 
sion must correspond to a luminosity near the Eddington limit, 
and we thereby confirm the distance to be close to 3 kpc. The av- 
erage accretion rate suggested by the broad-band non-burst flux 
is 0.32-1.6% of the Eddington rate. 

The two intermediate-duration bursts occur in a two-year 
time interval when the accretion flux lacks the otherwise so typ- 
ical flaring behaviour and is on average about 30% lower than at 
other times. Similar eff'ects are seen in other (candidate) UCXBs. 
A calm component is present which changes on time scales of a 
year, while flaring occurs with a time scale of a week (see in 't 
Zand et al. 2007). Also, bursts are shorter when a source shows 
more flaring activity (e.g., A 1246-58; in 't Zand et al. 2008). 
Moreover, UCXBs accreting above 1% of the Eddington limit 
show burst recurrence times on the order of days; those accret- 
ing below 1 % of the Eddington limit show recurrence times on 
the order of weeks (in 't Zand et al. 2007). In this respect it is in- 
teresting to note that during 4U 06 14+09 1 's high flux and flaring 
period it was detected in the radio (0.34±0.02 mJy at 4.86 GHz), 
while it was not detected during the low flux and calm period 
(3cr upper limit of about 0. 1 mJy at 4.86 GHz). This supports the 
existence of different states and suggests a connection between 
the radio and X-ray emission mechanisms (Migliari et al. 2009). 

The wide range of peak fluxes is not unprecedented. 
Galloway et al. (2008) find that roughly 2% of bursts have bolo- 
metric peak fluxes less than 0.1 times the maximum for the 
same source. Three sources exhibit dynamic ranges in excess of 
100: EXO0748-676 (236±36), 4U 1608-522 (446+115), and 
4U 1636-536 (151+25), although it should be mentioned that 
the high value in the last source is due to a weak secondary burst 
in a double burst (excluding that diminishes the dynamic range 
to 23+1). These three sources are not UCXBs, implying that the 
wide dynamic range in 4U 0614+091 is not necessarily due to 
extraordinary abundances of H, He, C or O. 

We found two long-lasting faint tails wherein the flux de- 
cays exponentially with a time constant in the range of 1000 
to 2000 s. One of them is the remainder of the intermediate- 
duration radius-expansion burst. The tail is seen clearly up to 
at least 2.5 hours after the start of the burst, possibly even up to 
5 hours. During the tail no cooling is seen; the inferred temper- 
ature stays constant near about 1 keV, while the apparent emit- 
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ting area decreases. Similar faint tails have been seen in various 
other bursters (Linares et al. 2009, Falanga et al. 2009, in 't Zand 
et al. 2009), but they were not as long lasting as those seen in 
4U 06 14+091. A hot underlying neutron star (with temperatures 
between 0.5-1 keV, i.e., slightly lower than that found in the 
burst tail) can explain the fact that in the tails the inferred tem- 
peratures are effectively constant (van Paradijs & Lewin 1986, 
in 't Zand et al. 2009). In this regard, we note that a soft compo- 
nent in our persistent PCA spectral fits is present, perhaps com- 
ing from the neutron star and/or the inner parts of the accretion 
disk; when the soft component is modelled by a black-body we 
infer temperatures of about 1 keV, i.e., close to that found in the 
long-lasting tails. As noted by in 't Zand et al. (2009), neutron 
stars in UCXBs accrete at lower rates and lack the CNO burning 
of H in their surface layers, compared to those in other, ordi- 
nary, bursters. This may result in a cooler neutron star in the for- 
mer systems, leading to a different behaviour in the long-lasting 
tails during their bursts. However, the temperature reached in the 
long-lasting faint tail of 4U 06 14+091 (and that of the transient 
XTE J 170 1-407, see Linares et al. 2009, Falanga et al. 2009, 
presumably also containing a cool neutron star due to its long 
off-states) is similar to that seen in ordinary persistent bursters, 
which may indicate that a different explanation is required for 
the constant temperature in the tail. 

4.2. The burst light curves arid energetics 

In an attempt to constrain the energy release and ignition 
depth of the bursts, we made some models of the light curves 
following the approach of Gumming & Macbeth (2004) for 
superbursts. We deposit an energy per unit mass E-En x 
lO'^ergg"' uniformly in the layer down to a particular col- 
umn depth y. The cooling of the layer is then followed by in- 
tegrating the thermal diffusion equation. We set the compo- 
sition after burning to be ^''Fe. Our models assume neutron 
star parameters of 1.4 Mq, /?=10km, giving a redshift factor 
l+z=1.31 and ^=2.44x10'^* cm s"^. To avoid following the ex- 
pansion of the outer layers that occurs when the luminosity 
approaches the Eddington limit, we set the top of our grid to 
y^lO'^gcm"^, and choose our outer boundary condition such 
that the maximum flux at that point is the Eddington flux at 
the surface of the star for pure He composition, i.e., F-Edd-cglK 
with /(:=0.2cm^g ' (Thomson scattering for pure He), giv- 
ing fEdd,oo=/^Edd/(l+z)^=2.1xlO^^ ergcm"^ s"'. Because of this 
crude treatment of the outer layers, and the simplistic way in 
which the energy is initially deposited, our models are not valid 
for early times, <10s (the thermal time scale at a depth of 
lO^gcm"^). 

We varied the two parameters £17 and y to find values that 
agree well with the observed light curves. Note that the values of 
£17 and y are taken to be independent parameters of the model, 
which allows us to investigate the ignition depth and energy re- 
lease without assuming a particular fuel for the burst. The model 
light curves are compared with the observed light curves for the 
superburst, the intermediate-duration burst on MJD51944, and 
the short burst on MJD 53718, in Fig. [10] We assume a distance 
to the source of d-3 kpc and a neutron star radius /?=10km to 
convert the observed fluxes into flux at the surface of the neutron 
star; these values give good agreement with the maximum flux 
from the model. 

The best fitting column depths for the three bursts in or- 
der of decreasing duration are y^ 2xl0"gcm-2, :^8xlO'' gem , 
and ^2x10^ gcm"^. For the intermediate-duration and short 
bursts, the column depths inferred from the light curves are 
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Fig. 10. Models of the burst light curves for the superburst (SB), 
the intermediate-duration burst on MJD 51944 (H1488), and 
flie short burst on MJD 53718 (H3982). For the superburst, the 
model parameters are £17=5, y\\-2 (solid curve) and En-2, 
3'ii=3 (dotted curve). For the intermediate-duration burst, the 
model parameters are En-l, ^9=8 (solid curve) and £17=5, 
3^9=10 (dotted curve). For the short burst, the model param- 
eters are £17=6, ^9=2 (solid curve) and £17=4, ^9=3 (dotted 
curve). Here E-Eii x lO'^ergg"' andy-ytj x 10*^ g cm"^=3'n x 
10" g cm"^. The observed fluxes and model fluxes for the short 
burst have been divided by a factor of ten for clarity. The start 
time of the superburst is unknown, but we find that choosing the 
start time to be around 4 hr before the time of the first data point 
gives a slope for the observed light curve that agrees with the 
models. Here we adopt a start time of 3.9 hr before the first data 
point. 



consistent with accretion at about 1% Eddington for the ob- 
served times since the previous burst, Af=155 and 16 days, 
respectively (see Table [T]i. This can be derived as follows. 
The local Eddington accretion rate (accretion rate per unit 
area) at the surface of the neutron star is riiEdd-c/RK. We 
define the Eddington rate as the value corresponding to 
/?=10km and /(■=0.2 cm^ g"', giving mEdd=l-5xlO^ g cm"^ s"', 
or 1% mEdd=1500g cm"^ s"'. Accretion at m=l%mEdd then 
gives a column 3'=mAf/(l+z)= 1.6x10'' gcm"^ (Af/16 d)(m/l% 
'WEdd)- Matching the column depths infeiTed from the light 
curves assuming that At is the recurrence time implies m=0.52%, 
and 1.3% mEdd for the intermediate-duration and short burst, re- 
spectively. 

The time span between the superburst and the previous ob- 
served burst is 367 days (i.e., Af=367 days, see Table [U, in 
which time the accumulated column is y = ihAt/{l + z) = 
3.7 X 10"' g cm"2 (Af/367 d)(m/l% mEdd)- Therefore, we require 
the accretion rate to be m^5% mEdd for the column infeiTed from 
the superburst light curve to be accreted in the 367 days leading 
up to the superburst (assuming no bursts occurred in between). 
An alternative explanation is that the superburst involves a dif- 
ferent fuel layer (e.g.. He for the short and intermediate-duration 
bursts, C for the superburst) in which case the column for the 
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Fig. 11. The observed burst energies (assuming c/=3kpc) ver- 
sus the time since the last observed burst Af from Table [U 
Because of possible missed bursts due to the many data gaps, 
Af is an upper limit on the average burst recurrence time. The 
solid lines show the expected relation for accretion at 1% of 
the Eddington rate assuming complete burning of the accumu- 
lated fuel for two different values of nuclear energy release per 
nucleon, 2,iu(; = 1.6MeV coiTesponding to complete burning of 
pure He to Fe group, and a smaller value 2nuc=0.6 MeV. The 
dotted, dashed, and dot-dashed curves show the predicted rela- 
tion for pure He accretion from Gumming et al. (2006). (The 
different curves are for different core neutrino emissivities and 
crust properties; see Gumming et al. [2006] for details.) 



superburst could accumulate over a longer time scale than Af. At 
1 % Eddington, the accumulation time would then be about 5 yr. 

The best fitting values for the energy release are En^S-l 
(see Fig. [TOb . In addition, the fact that the light curves reach 
the Eddington flux implies a minimum energy release. For the 
intermediate-duration burst, we find that the flux does not reach 
the Eddington flux for En<4. These values of energy release are 
smaller than the values £'i7a;16 for complete burning of He to 
Fe group or En^lO for complete burning of G to Fe group. This 
implies that the energy release is about a factor of 3 lower than 
expected for complete burning of pure He. 

The low energies of the bursts can also be seen in the over- 
all energetics. If all the matter accreted since the previous burst 
were completely burned during each burst, one would expect the 
ratio, a, of the average luminosity in the persistent emission to 
the time-averaged luminosity emitted in type-I X-ray bursts to 
have a value in the range 25-200, depending on the composi- 
tion of the burning material. If one assumes that both the persis- 
tent and burst emission are isotropic, that the persistent emis- 
sion has not varied since the previous burst, and that Fp^-s is 
the bolometric persistent flux, then this statement relating the 
persistent and burst luminosities is equivalent to Q;=AfFpers/£b- 
If we take that the minimum average burst recurrence time is 
the lowest observed value for Af, i.e., 7 days, then the corre- 
sponding values for a are in fact lower limits, which are be- 
tween 125 and 5525 for the short bursts and around 30-70 for 



the intermediate-duration bursts. Galculating an a value by as- 
suming that no bursts were missed between observed bursts (that 
is Af is the average burst recurrence time) gives a range from 
500 to 5000, larger than the canonical value of 100-200 for pure 
He burning. Another way to look at this is presented in Fig. [TT] 
which shows the burst energy against time since the last burst 
Af from Table [1] The solid lines show the expected burst en- 
ergy as a function of recurrence time for accretion at 1 % of the 
Eddington rate, for two different values of the nuclear energy 
release per nucleon 2nuc=l-6 and 0.6MeV per nucleon (where 
£'i7~10(2„uc/MeVpernucleon)). The curves at the top left are 
pure He ignition models from Gumming et al. (2006) that as- 
sume complete burning. The observed bursts lie well below the 
expected energies assuming complete burning. 

Fig. [To] shows also that for the intermediate-duration burst 
on MJD51944, the model light curve nicely connects the early 
decay from the peak luminosity to the long tail lasting thousands 
of seconds. Therefore the long tail is naturally explained by the 
cooling of deep layers heated by the burst as proposed for tails 
seen in other sources by in 't Zand et al. (2009). The theoretical 
expectation is that the luminosity should decay as a power law 
in the cooling tail. This is confirmed by our observations (e.g., 
Fig.Iinil. 

4.3. He as a fuel for the bursts 

In Sect. 14.21 we found that the short and intermediate-duration 
burst light curves could be explained by ignition depths of y^2 
and 8xlO''gcm"^, respectively, and that these depths are com- 
patible with the amount of mass accreted at an accretion rate 
close to ih^l% riiEdd in the observed time between bursts, Af. 
These ignition depths would be naturally explained if the neu- 
tron star accreted a substantial amount of He. At low accretion 
rates, most of the energy released by pycnonuclear and elec- 
tron capture reactions in the crust flows outwards, because the 
neutron star core is too cold for significant neutrino emission. 
This results in an outward flux from the crust of F-ihQh with 
2i!,^lMeV per nucleon (Brown 2000; Gumming et al. 2006: 
their figure 18, top panel). For this value of Qh, figure 22 of 
Gumming et al. (2006) shows that pure He accretion gives igni- 
tion aty~10'^-10"'gcm"^ at m~l%m^i. 

However, as we noted in the introduction, the lack of He 
lines in the optical spectrum of the source implies that the ac- 
creted material is not significantly enriched in He. To investigate 
whether a small amount of He in the accreted material could still 
explain the observed bursts, we constructed ignition models for 
type I X-ray bursts following Gumming & Bildsten (2000). We 
assume that the accreted material consists of He and a 50/50 G-O 
mixture. We find that for small He mass fractions, the He burns 
away stably before reaching the conditions required for a ther- 
mal runaway. In Fig.[T2]we plot the He mass fraction required to 
achieve unstable ignition as a function of m and Qi,. It shows that 
the amount of He inferred from the fits of disk spectra by Werner 
et al. (2006), i.e., y<10%, do not lead to unstable ignition at any 
accretion rate unless the neutron star is cold, Qh^Q. 1 MeV per 
nucleon, and even then the recurrence times are extremely long. 
We, therefore, conclude that matching the observations requires 
a significant amount of He in the accreted material. 

If we allow a large enough He fraction, we can match 
the inferred ignition depth for the short and intermediate- 
duration bursts. For example, taking a He mass fraction 7=0.5, 
2fc=0.3 MeV per nucleon, and local accretion rates m=3000 and 
5000gcm"^ s"' (2% and 3.3% of the Eddington rate) gives ig- 
nition depths of y-2 and 8x10^ g cm"^, respectively. The coiTe- 
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Fig. 12. The mass fraction of He needed to achieve unstable He 
ignition as a function of accretion rate, for accretion of He and an 
equal mixture of C and O. The different curves are for different 
values of Qb as labeled. 



sponding recurrence times are 7 days and 46 days, respectively. 
For pure He and Qh-iyitW per nucleon, the correct ignition 
depths are obtained for m-\% and 1.7%mEdd, with recurrence 
times of 81 days and 11 days for the intermediate-duration and 
short burst, respectively. These recuiTence time estimates are 
consistent with the observed values. 

For each of these examples, a decrease in accretion rate of 
only 40% is required to change the ignition column depth by a 
factor of 4 from 2 to 8x10^ g cm"^. This is because of the steep 
dependence of the ignition column on temperature - only a small 
decrease in temperature, and a small corresponding decrease in 
m is needed. The required decrease in accretion rate is compa- 
rable to the observed difference in persistent flux between the 
2001-2002 calm period when the intermediate-duration bursts 
occurred and the 2003-2007 flaring period when more regu- 
lar bursting resumed: as noted in Sect. 14.11 the 1 -year-averaged 
ASM fluxes for 4U 0614+091 in 2001 and 2002 are about 30% 
smaller than for the years after this time interval. 

It is not clear why unstable He ignition would give a nu- 
clear energy release of only 2nuc~0.6 MeV, as inferred from the 
burst energetics and light curve fits (see Sect. I4.2l i. This en- 
ergy is approximately the energy released in burning He to C. 
For a He mass fraction Y, setting the total energy release to be 
Gnuc=0.6 MeV per nucleon implies that the energy release from 
burning beyond C to heavy elements is 2heavy~0.6(l - y)MeV 
per nucleon. This is much less than the » 1 MeV available for 
complete burning of C to Fe, unless the helium fraction Y is 
small (however, as noted above, the observations require Y to be 
large). At a depth of y~\(f g cmT^, an energy deposition of only 
0.03 MeV per nucleon is sufficient to raise the temperature to 
above 10^ K, at which the nuclear burning would be expected to 
proceed beyond C, so a significant energy release from burning 
beyond C would be expected. Numerical models of bursts at low 
accretion rates would be useful to follow the nucleosynthesis and 
determine the expected energy release. 



Another possibility is that not all of the accreted fuel burns 
during the burst. For example, the burning may consume only 
part of the depth of the fuel layer, or may cover only part of the 
stellar surface. Numerical models of bursts at low accretion rates 
are needed to follow the nucleosynthesis and burning dynamics, 
and determine the expected energy release. 

4.4. The superburst 

In the previous section we argued that the short bursts and 
intermediate-duration bursts could be explained by a substantial 
fraction of He in the accreted material, despite the fact that the 
optical spectrum of 4U 0614+091 suggests that very little He is 
present in the accreted material. However, the superburst poses 
a more severe problem. 

The superbursts observed in other sources at accretion 
rates m>Q. 1 medd have been explained as being due to C igni- 
tion (Strohmayer & Brown 2002, Gumming & Bildsten 2001). 
However, there are problems with this scenario: (1) producing 
enough C during H/He burning (e.g., Schatz et al. 2003, Woosley 
et al. 2004), (2) heating the neutron star ocean strongly enough 
to reach ignition temperature (e.g.. Gumming et al. 2006, Keek 
et al. 2008), and (3) accreting rapidly enough for the C to survive 
to the ignition depth (Gumming & Bildsten 2001, Gumming et 
al. 2006). 

If there is accretion of a significant amount of G from the 
GO white dwarf companion in 4U 0614+091, the first problem 
would be eased by removing the need to make the G during nu- 
clear burning of H or He. However, if He flashes are responsi- 
ble for the intermediate-duration and short bursts, the G would 
likely burn away during these flashes. The last two problems are 
more difficult to circumvent. First, the ignition temperature for 
G at the column depth inferred from the superburst light curve, 
37^2x10" g cm"^, is above 6x10^ K. This temperature is much 
higher than the ignition temperatures of He flashes, i.e., about 
1x10*^ K. Achieving this temperature is even more difficult when 
the accretion rate and therefore crust heating rate are signifi- 
cantly below those typical around the times of the superbursts 
seen in other sources. Second, we calculated models for G igni- 
tion following Gumming & Bildsten (2001) and extended their 
figure 2 to lower accretion rates, but were not able to find values 
of Qh or G fraction for which the G survives to ignition depth at 
accretion rates of ~0.01 mEdd- 

Another possibility is that ignition of a thick He layer 
is responsible for the observed superburst. Indeed, for the 
Qb=0.3 MeV per nucleon and 7=0.5 case considered earlier, we 
find that reducing the local accretion rate to 1000 g cm"^ s"', a 
factor of 3 below the accretion rate that reproduces the igni- 
tion column of the intermediate-duration burst, gives He ignition 
at 2x10" gcm"^. The time scale for He burning is longer than 
the accumulation time (although close to it), indicating that He 
should survive down to the ignition depth, even for such a low 
accretion rate (this is not the case for the Y-l, Qb-l MeV per 
nucleon model; there we find that the He burns stably away). 
However, although the ignition depth can be achieved, the ex- 
pected recurrence time is 9.8 yr, an order of magnitude longer 
than the observed time between the superburst and the previous 
burst, Af=367 days. On the other hand, this can be reconciled if 
an appreciable portion of the accreted He survives the preceding 
shorter bursts (see previous Section). 

Another constraint on the superburst ignition depth comes 
from the quenching time scale for normal bursts following the 
superburst. Normal bursting behaviour resumed 19 days follow- 
ing the superburst, the fastest time scale so far observed (see. 
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e.g., Kuulkers 2004). Gumming & Macbeth (2004) predict that 
the quench time for normal bursting behaviour after a super- 
burst is fquench=l-6yi2'*(OT/OTEdd)"^^'*£'if days (where y^yu X 
10'^ g cm"^). The constraint on the quench time rquench<19 days 
implies for En -6 and (m/mEdd)=0.016 (at the time of the su- 
perburst) that >'i2<0.18. This is in good agreement with the ig- 
nition column inferred from the superburst light curve, and con- 
sistent with the normal burst observed 19 days after the super- 
burst being the first burst to occur following the superburst. The 
dependencies are not very strong. For M/MEdd=0. 0032-0.016, 
the ignition depth is between )'i2<0.04-0.18. For £17=1-6 this 
becomes >'i2<0. 18-0.44. Therefore, the quench time constraint 
provides another piece of evidence that the ignition depth for 
the superburst was lower than the depths inferred for all or most 
previously analysed superbursts in other sources. 

5. Conclusions 

We can understand several aspects of the type I X-ray bursts ob- 
served from 4U0614-H091. The colurmi depths and energy re- 
lease per gram in the bursts can be estimated by comparing the 
observed light curves with models, and by considering the burst 
energetics. Both methods are in good agreement. Helium igni- 
tion naturally explains the observed ignition depths for accretion 
at m~l% ffiEdd- Furthermore, the sensitive dependence of the He 
ignition depth on temperature means that small (factor of two) 
changes in accretion rate can lead to an order of magnitude vari- 
ation among the depths, and can thereby explain the occurrence 
of both short bursts and intermediate-duration bursts. The igni- 
tion depth for the superburst inferred from the Ught curve is the 
lowest of the current sample of superbursts and agrees well with 
the constraint from the observed quenching time scale. 

However, several puzzles remain. First, the amount of He re- 
quired to achieve the required ignition conditions without sta- 
bly burning away is significantly larger than the ^10% limit 
from the optical spectra. Recently, however, formation stud- 
ies using evolutionary calculations indicate that the donor in 
4U 0614-H091 may be a hybrid white dwarf or very evolved he- 
lium star (Nelemans et al. 2009). This suggests the donor still to 
be a possible supplier of a significant amount of He to the accre- 
tion disk, onto the neutron star. Further investigations have to be 
done, why the He does not show up then in the optical observa- 
tions. 

Second, understanding the superburst remains problematic. 
Unstable C ignition is not known to be possible at these low 
accretion rates: the accumulating layer is too cold, and even if 
heated the C burns stably. A superburst powered by a large He 
pile is possible, but takes several years to accumulate if the ac- 
cretion rate is ~1% Eddington. Such an accumulation time is 
much greater than the observed time of one year Finally, the 
reason for the low energy per gram ^nuc~0.6MeV per nucleon 
released in the bursts is not clear. 

It is important to emphasize, however, that the ignition mod- 
els and light curve models used in this paper are simplified. The 
model light curves assume uniform and instantaneous energy de- 
position in the fuel layer, and do not follow the detailed nucle- 
osynthesis. These models cannot address the early part of the 
light curve, for example the interesting ripples in the early light 
curve of the intermediate duration burst (Fig. 10), nor whether 
the burning is expected to be incomplete (as we infer from the 
burst energetics). Incomplete burning, either across the neutron 
star surface or with depth in the fuel layer, may at least partly 
solve the above puzzles. Our ignition models rely on a simphfied 



one-zone ignition criterion. Numerical models of accumulating 
and burning of He/CO mixtures at low accretion rates are needed 
to confirm our conclusions, for example, regarding the amount 
of helium needed to avoid stable burning. 

It is interesting that at least in principle He can power a 
superburst-like event. The requirement is that the accumulating 
fuel layer remains cold. Recently, Cooper et al. (2009) ruled out 
He as a fuel for superbursts, but their argument assumed an igni- 
tion temperature of 7=5x10^ K that is much larger than the He 
ignition temperature at superburst columns. Even in the super- 
burst sources with m>0. 1 mEdd, He -powered superbursts could 
occur if the accumulating layer is cold enough. One way to keep 
the layer cold would be to have direct URCA neutrino emission 
in the neutron star core, so that most of the energy release in 
the crust flows inwards rather than outwards. Triple alpha igni- 
tion becomes mostly sensitive to density rather than temperature 
when the colunm depth reaches y~3y.\Q^^ g cm"^, as the ignition 
becomes pycnonuclear and therefore temperature independent. 
This could potentially explain the narrow range of superburst ig- 
nition columns. Further work on this possibility is needed. 
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Appendix A: Localizations of bursts 

The three bursts observed by EURECAj^NlCW had positional 
3cr error circles of less than 1 ° radius (taking into account sys- 
tematic errors and errors due to uncertainties in the spacecraft 
pointing). The optical (Davidsen et al. 1974, Murdin et al. 1974), 
infrared (Migliari et al. 2006) and radio (Migliari et al. 2009) 
counterparts to 4U 0614-1-091 were always within these error cir- 
cles. The WATCH imaging results for the event which occurred 
on February 17 is shown in Fig |A.ll (left)F^ 

At the time of the sole burst seen by the WFC, the satellite 
attitude solution was not optimum. However, the position for the 
burst is the same to within 0. 1 pixel of the position determined 
for the persistent emission at times when the attitude solution 
was optimum, with an uncertainty of 2' (99% confidence). We 
conclude that the origin of the burst is coincident within 2' of 
the position of 4U 06 14-H09 1 . 

The first two bursts (in 1996 and 1998) of the six normal 
ASM bursts were observed by multiple SSCs, either because the 

Note that this analysis was done more than 10 years ago. Therefore, 
the best-fit estimates on the position of the various bursts are not avail- 
able; also the regeneration of images is not possible anymore. 
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Fig.A.l. Left: The cross-correlation map (6-15 keV) contain- 
ing the burst observed from the direction of 4U 0614-1-091 on 
February 17, 1993 by EURECAI^NPC^CW. The burst was con- 
tained in a modulation pattern integrated for 56.7 s. The con- 
tours are shown in Crab units and the averaged correlation at 
the position of 4U 0614+091 was 3.6 Crab. Note that the in- 
tegration time is too short to give a detection of the Crab. 
Right: Detection of the persistent emission from 4U 06 14+091 
by Et/RZiCA/WATCH. The skymap is based on data from the 
period of February to April 1993 in the energy band 6-12 keV. 
The contours show the cross correlation in units of mCrab (only 
positive coiTelation is shown). The ring-like structure around 
4U 0614+091 is an artefact caused by the WATCH RMC detec- 
tion principle. 

































'I>>1^ 




























■ \ \ \ s 


^ / 


fiiiOQ" 






.^^^ / 











Fig. A.2. Equatorial (J2000.0) map of ten independent ASM lo- 
calizations of the source of six bursts. The asterisk marks the 
location of the optical counterpart to 4U 06 14+091 (Davidson et 
al. 1974, Murdin et al. 1974). 



burst position was in the region where the fields of view of SSC 
and SSC 1 overlap, or because the burst was active across two 
sequential dwells (or both, as in the first case, yielding four in- 
dependent observations). The four subsequent bursts were only 
detected in a single SSC for a single dwell, yielding a total of 
eleven observations of the six bursts (see TableO. 

For each ASM dwell, the intensities of known sources in the 
field of view are derived via a fit of model slit-mask shadow 
patterns to counts binned by position along each anode in each 
detector. The residuals from a successful fit are then cross- 
correlated with each of the expected shadow patterns corre- 
sponding to a set of possible source directions which make up 
a grid covering the field of view. A peak in the resulting cross- 



Table A.l. Observation log of ASM bursts from 4U 0614+091 
(left) and 2S 0918-549 (right).^ 



4U 0614+091 2S 0918-549 



MJD 


SSC 


fx 


MJD 


SSC 


fx 


50200 





0.648/0.419 


50211 





0.418 




1 


0.536/0.522 


52181 


1 


0.661 


51164 


0,1 


0.648,0.147 


52509 





0.208 


53476 





0.512 


52856 


0,1 


0.698, 0.230 


53959 





0.588 


53683 


1 


0.671 


54101 


2 


0.141/0.044 








54108 


1 


0.142 









Note 3. Given are the day (MJD) at which the burst occurred, with 
which SSC it was seen (0 or 1) and the transmission factors (fx) for 
the corresponding SSCs (X=0 or 1) in the 1.5-12keV band. For the 
first and fifth bursts of 4U 0614+091 we give the transmission factors 
for two sequential dwells, separated by a slash. 



Table A.2. Positions of the bursts seen by ISGRI and BAT. 



Instrument/ 


RA 


Dec 


error" 


MJD 


(J2000.0) 




ISGRI 








52867 


94.281° 


9.152° 


2.1' 


53460 


94.284° 


9.133° 


7.0' 


BAT 








54029 


94.278° 


9.149° 


1.5' 


54189 


94.203° 


9.155° 


3.1' 



90% confidence error region. 



coiTelation map indicates the possible presence and approximate 
location of a new, uncatalogued X-ray source (Levine et al. 1996, 
Smith et al. 1999). 

To test the hypothesis that the six ASM bursts were from 
4U 06 14+091, we removed this source from the catalog used for 
the fitting procedure described. We fit the eleven observations 
with this truncated catalog, and we searched the residuals for 
evidence of a source of emission unaccounted for by the cata- 
log. For bursts that occurred after the ASM began recording in 
event mode, we omitted from the fit those time bins that did not 
show evidence of burst activity, to increase the signal to noise 
ratio. All eleven dwells showed evidence of a 'new' source in 
the residuals. If 4U 06 14+091 were the source of the burst in 
each of these dwells, the localization of this new source should 
be consistent with the location of 4U 0614+091. 

Ten of the eleven localizations are plotted in Fig. IA.2I The 
boxes are determined by the 95% confidence interval in each 
of the two dimensions recorded by the SSC, projected onto the 
sky. Since the two dimensions are independent, the total confi- 
dence level of the eiTor box is about 90%. These uncertainties 
are statistical only. The location of the optical counterpart of 
4U 06 14+091 is indicated by an asterisk. This location cannot 
be excluded by any of the error boxes, including an eleventh er- 
ror box that lies outside the limits of this plot. 

For ten of the thirteen FREGATE bursts a precise localiza- 
tion was possible using the WXM and/or SXC, and they are con- 
sistent with coming from 4U 0614+091 (see Sect. 2.1.5). For the 
other three bursts we can only suppose they are coming from 
4U 06 14+09 1 , since the source was the only persistent burster in 
the FREGATE field of view. 

The reconstructed 15-30keV images of the field of view 
around 4U 06 14+09 1 during the bursts seen by ISGRI are shown 
in Fig. IA.3I The only two detected sources in the 2003 (24 s) 
and 2005 (15 s) time frames are indeed 4U 0614+091, as well 
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Fig.A.3. IBIS/ISGRI (15-30keV) significance images during 
the bursts detected on 2003 August 16 {top) and 2005 March 31 
(bottom). Shown is the equatorial J2000.0 grid with 5° spacing. 
4U 06 14+091 and Crab are the only sources detected in the time 
frames (UT 21:30:18 - 21:30:43 and UT 07:12:18 - 07:12:34, 
respectively). The detection significances of 4U 06 14+091 (and 
Crab) in the 2003 and 2005 time frames ai'e 4.7 (18.3) and 12.7 
(6.9), respectively. 



as the persistently bright Crab source. The data from the two 
BAT bursts were subjected to the post 'y-ray burst' processing 
script, which led to refined positions for these bursts. The derived 
ISGRI and BAT bursts coordinates are given in Table I A. 2] They 
are also consistent with the optical, IR and radio counterparts to 
4U 06 14+091. 

We conclude that the 27 bursts seen by Ef/RECA/WATCH, 
WFC, ASM, FREGATE, ISGRI and BAT indeed originated from 
4U 06 14+091. 



